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EXECUTIVE  PROGRAM  SUMMARY 


A  major  accomplishment  of  this  year's  NRL  Plasma  Radiation  Source  (PRS) 
Theory  Program  was  a  demonstration  of  how  long-term  theoretical  support  allows 
sparse  and  diverse  experimental  results  to  be  organized,  interpreted,  and 
understood.  As  a  result,  guidance  that  will  lead  to  improved  simulation 
fidelity  was  given  this  year  for  the  design  of  pulsed  power  machines  as  plasma 
radiation  sources  (PRS).  To  accomplish  this,  a  large  series  of  pre-planned 
wire  array  aluminum  hydrodynamics  calculations  were  carried  out  and  the 
results  were  extensively  analyzed.  This  analysis  provided  a  set  of  minimum 
requirements  for  pulse  power  machine  design.  It  predicted  aluminum  K-shell 
yield  behavior  as  a  function  of  array  mass  and  current,  led  to  predictions  of 
how  a  given  K-shell  yield  scales  with  atomic  number,  and  stimulated  a  new 
series  of  PRS  experiments  that  were  designed  and  carried  out  at  Physics 
International  Inc.  using  aluminum  wire  arrays.  These  experiments  will  be 
analyzed  in  the  coming  year.  Preliminary  analysis  indicates  that  the 
calculations  had  indeed  predicted  minimum  aluminum  K-shell  yields,  and  it  led 
to  the  observation  that  at  least  one  other  important  x-ray  energy  conversion 
mechanism  is  at  work  in  the  PRS  plasma.  One  such  mechanism  has  received  some 
preliminary  phenomenological  attention  and  will  be  described  later  in  this 
final  report.  A  data  base  needed  to  analyze  PRS  aluminum  plasmas  for  their 
temperature,  density,  and  K-shell  mass  participation  was  also  completed.  It 
will  allow  a  two  pronged  analysis  of  the  PI  PRS  data  to  be  carried  out  and 
quantities  such  as  anomalous  resistivity  and  other  anomalous  plasma  properties 
to  be  unambiguously  inferred  from  the  data. 

Other  important  NRL  theoretical  efforts  continued  this  year.  They  fall 
into  five  program  categories: 

I.  Program  Technical  Guidance  and  Data  Analysis 

o  Detailed  K-shell  emission  atomic  number  scaling  laws  were  developed 
for  kinetic  energy  to  x-ray  conversion  and  are  being  evaluated  using  nickel 
array  implosion  calculations. 


o  Maxwell  Laboratory  aluminum  array  implosion  data  were  analyzed  and 
trends  in  plasma  temperature,  density,  and  mass  participation  were  determined. 

o  The  approach  to  equilibrium  in  a  reflex  switch  was  analyzed.  The 
predicted  and  the  observed  voltage  and  current  waveforms  were  found  to  be  in 
good  agreement  so  there  is  hope  that  some  basis  for  understanding  reflex 
switch  operation  has  been  established. 

II .  Atomic  Data  Acquisition 

o  A  complete  set  of  atomic  data  was  calculated  for  F  to  Ne  dielectronic 
recombinations  in  four  different  elements.  This  data  was  found  to  be 
scalable.  Consequently  it  will  be  used  to  produce  reliable  DR  data  for  Ar 
through  Kr  ions,  and  it  will  lead  to  reliable  satellite  line  spectral 
diagnostics . 

o  A  complete  collision  strength  data  base  for  calculating  Fe  and  Se  L- 
shell  collisional  excitations  was  calculated.  Preliminary  analysis  suggests 
that  this  data  should  be  scalable.  This  means  we  will  get  20  data  sets  for 
the  price  of  3. 

o  Procedures  for  calculating  and  scaling  radiative  decay  rates  in  the 
L-shell  were  also  developed. 

III .  Radiative  Collapse 

o  Conventional  MHD  descriptions  of  an  imploding  alumint  a  (or  low  Z) 
array  using  present-day  pulse  power  voltage  waveforms  generally  lead  to  a 
radiative  collapse  in  contradistinction  to  the  experimental  situation.  A 
careful,  phenomenological  analysis  of  PI  aluminum  experiments  was  carried  out 
to  demonstrate  that  a  short  circuit  in  the  diode  and  anomalous  resistivity  in 
the  plasma  will  move  the  calculated  x-ray  emission  in  the  direction  of 
experiment  and  eventually  produce  fair  agreement. 


v 


IV. 


PRS  Optimization 


o  The  problem  of  determining  machine  coupling  efficiency  to  future  PRS 
loads  was  addressed  in  three  separate  ways: 

o  The  way  a  diode  short  circuit  affects  machine  coupling  in  an  argon 
puff  gas  implosion  is  currently  being  investigated. 

o  The  bleeding  of  current  away  from  the  main  body  of  plasma  by  a  plasma 
corona  was  investigated,  and  the  theory  of  this  corona  was  critiqued.  In 
general,  we  need  to  evaluate  the  degree  of  machine  PRS  decoupling  produced  by 
runaway  electrons,  and  the  existence  of  a  corona  may  be  a  relevant  aspect  of 
this  problem,  since,  it  was  determined  this  year  that,  our  present-day  MHD  is 
not  in  accord  with  existing  experimental  data  on  runaway  electron  production. 

o  The  coupling  efficiency  of  a  transmission  line  to  screwpinch  and  z- 
pinch  PRS  loads  was  calculated  and  found  to  be  acceptably  large. 

V.  New  PRS  Approaches 

o  By  imploding  onto  and  compressing  a  magnetic  field  in  a  screwpinch 
array  implosion,  one  can,  in  principle,  control  implosion  stability  in  the 
laboratory.  An  ideal  MHD  screwpinch  model  was  developed  this  year  and 
benchmarked  successfully  against  self-similar  analytical  MHD  flows.  Once 
magnetized  transport  and  magnetic  field  diffusion  are  added  to  the  model,  it 
will  be  used  to  carry  out  detailed  calculations  of  the  stability  of  presently 
achievable  magnetic  shear  profiles. 

o  Ring  plasma  stagnation  dynamics  between  an  argon  plasma  and  an  iron 
wall  were  calculated.  They  confirmed  that  the  stagnation  of  a  high  velocity 
ring-puff  plasma  on  a  planar  solid  density  target  will  generate  a  high  density 
plasma  that  can  efficiently  emit  multi-keV  radiation. 


o  A  2-d  MHD  model  of  a  copper  x-pinch  was  developed  to  investigate  the 
dynamical  evolution  of  the  pinch  and  to  predict  the  variation  in  radiated 
energy  as  a  function  of  the  mass  per  unit  length.  This  pinch  is  of  special 
interest  because  it  produces  a  point  source  of  x-rays. 

A  more  detailed  description  of  the  above  work  can  be  found  in  Sections  I 
through  V  of  this  report,  which  correspond  directly  with  the  5  program 
categories  listed  above.  Further  details  have  already  been  given  in  a  number 
of  papers  that  were  written  this  year  for  publication. 

As  a  result  of  the  work  that  was,  and  is  being,  carried  out  in  the  NRL 
Theory  Program  this  year,  a  number  of  challenging  new  ideas  have  emerged  for 
(1)  making  machine  upgrades  to  increase  K-shell  output  on  existing  PRS  pulse- 
power  machines,  (2)  carrying  out  new  experiments  with  more  efficient  load 
designs,  (3)  improving  our  understanding  of  present-day  machines,  and  (4) 
exploring  the  possibility  of  achieving  radiative  collapses.  For  example,  it 
is  anticipated  that  a  proper  upgrade  of  the  Phoenix  PRS  machine  at  NSWC  should 
produce  kilojoules  of  nickel  K-shell  radiation.  Other  new  developments  in  the 
status  of  the  PRS  program  can  be  summarized  as  follows: 

o  We  now  have  the  capability  to  learn  a  great  deal  about  PRS  loads 
using  time  resolved  experimental  diagnostics.  We  also  have  the  tools  to 
analyze  these  experiments  and  clearly  delineate  anomalous  heating  and  mass 
participation  effects.  These  diagnostics  and  tools  need  to  be  exercised. 

o  A  systematic  new  theoretical  interpretation  can  now  be  given  to  a 
wide  range  of  experimental  results  on  array  implosions.  However,  a  large  data 
base  of  such  results  needs  to  be  obtained  to  learn  more  of  the  phenomenology 
of  anomalous  heating.  This  development  is  of  special  significance  since,  in 
the  past,  PRS  experiments  conducted  on  a  given  machine  by  varying  load 
parameters  and  materials  have  been  notoriously  difficult  to  interpret  and 
misleading  when  it  comes  to  optimizing  PRS  design. 


o  Because  it  is  important  to  build  a  data  base  of  experimental  results 
from  as  many  different  PRS  machines  as  is  possible,  it  is  important  that  the 
decade  machine  be  designed  to  have  a  PRS  research  capability,  preferably  with 
options  for  variable  diode  and  current  design  (and  not  just  variable  load 
options) . 

o  We  now  have  a  credible  and  growing  theoretical  capability  to  predict 
mass  and  current  scalings  for  imploding  different  atomic  number  PRS  array 
designs.  However,  it  must  be  tested,  benchmarked,  and  improved  both 
experimentally  and  theoretically.  A  PRS  capability  on  the  decade  machine 
would  aid  and  abet  this  process. 

o  Most  importantly,  the  scaling  work  that  was  done  this  year  has  made 
it  clear  that  pulse  power  machines  and/or  their  diodes  must  be  designed  for 
the  specific  yield  and  spectrum  (or  element)  requirements  desired.  If  done, 
the  AGT  simulation  fidelity  goals  of  DNA  can  be  met  rather  than  approximated. 

o  Finally,  a  sound  and  viable  theoretical  foundation  has  been  built  to 
treat  problems  in  radiative  hydrodynamics  and  plasma  radiation  source 
development.  It  needs  continued  long  term  support  to  be  brought  to  full 
fruition  in  order  to  (1)  help  set  the  parameters  for  new  machine  design,  (2) 
aid  in  the  development  of  new  PRS  concepts,  (3)  provide  the  tools  necessary  to 
carry  out  data  analysis  and  interpretation  and  (4)  promote  the  design  and 
evaluation  of  new  experiments  and  new  improved  loads. 


ADVANCED  CONCEPTS  THEORY  ANNUAL  REPORT  1990 


I.  Program  Technical  Guidance  and  Data  Analysis 

The  analysis  of  the  scaling  of  K-shell  emission  in  aluminum  with  mass 
and  kinetic-energy-per-ion  was  completed.  In  this  work  we  investigated  a 
special  physical  situation  in  which  the  conversion  of  implosion  kinetic  energy 
into  K-shell  emission  at  stagnation  could  be  studied  unambiguously.  Some  of 
this  work  was  communicated  to  the  z-pinch  community  at  the  Laguna  meeting 
(April  26,  2989)  and  at  the  Colorado  Springs  Simulation  Fidelity  Workshop 
(September  19  -  23,  1989).  A  much  more  detailed  description  of  this  work  was 
given  in  two  papers  that  were  submitted  for  publication  in  JAP  and  JQSRT. 

In  the  JQSRT  paper,  the  mechanism  by  which  the  imploded  plasma  makes 
4  2 

the  transition  from  I  to  I  yield  scaling  is  described.  The  JAP  paper 
provides  a  complete  description  of  how  K-shell  emission  in  aluminum  scales 
with  mass  and  kinetic-energy-per-ion,  and  it  also  makes  predictions  as  to  how 
these  results  should  scale  with  atomic  number.  We  found  that  under  optimum 
conditions  about  30%  or  more  of  the  total  kinetic  energy  stored  in  the 
aluminum  plasma  can  be  emitted  as  K-shell  radiation.  The  atomic  number,  Z, 
scaling  was  done  by  assuming  that  this  conversion  efficiency  is  constant  for 
all  Z  materials.  However,  if  this  efficiency  is  substantially  less  for  higher 
Z  materials,  then  our  scaled  results  are  overly  optimistic.  We  have  initiated 
studies  using  nickel  in  order  to  determine  more  accurately  how  the  conversion 
efficiency  scales  with  Z. 

Chris  Deeney  at  PI  took  the  initiative  to  set  up  experiments  with 
imploding  aluminum  wire  arrays  to  put  the  above  aluminum  calculations  to 
experimental  test.  This  provided  us  with  an  excellent  opportunity  to  not  only 
work  and  exchange  ideas  with  the  people  at  PI  but  also  to  extend  the  aluminum 
scaling  predictions,  in  which  the  current  was  turned  off  just  prior  to 
implosion,  to  the  experimental  measurements  of  K-shell  emission,  where  the 
current  persists  throughout  the  implosion.  Since  the  purpose  of  the 
experiments  was  to  compare  experiment  and  theory,  maximum  machine  energy 
coupling  had  to  be  insured  by  having  the  implosion  times  occur  near  the  peak 
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in  the  open  circuit  current  profile  of  Double  Eagle.  According  to  0  and  1-D 

theory,  the  implosion  time  can  be  held  fixed  if  m  and  r  are  chosen  such  that 

2  2  ° 

(mrQ  )  *  x  ,  where  x  is  a  fixed  implosion  time.  Then  by  varying  m  and  rQ 

along  the  path  defined  by  mrQ2  =  constant  x  t2p,  where  is  the  implosion 

time  that  corresponds  to  the  peak  in  the  short  circuit  current  profile,  one 

can  determine  how  the  K-shell  emission  varies  with  mass  and  radius. 

Because  these  experiments  were  conducted  on  a  machine  with  a  given 

current  risetime  and  peak  current,  they  could  not  be  conducted  to  check 

against  the  calculations  directly,  since  the  calculations  were  carried  out 

holding  either  the  load  mass  fixed  or  the  kinetic-energy-per-ion  fixed. 

However,  using  the  current  peak  and  risetime  of  the  experiment  and  the 

condition  mrQ  =  constant,  one  could  theoretically  determine  the  mass, 

kinetic-energy-per-ion,  current,  and  implosion  time  variations  that  should 

have  occurred  in  the  experiments.  Then,  using  the  analytic  fits  to  the  MHD 

calculations  that  were  obtained  in  the  analysis,  one  could  determine  the  K- 

shell  yield  that  should  have  occurred  if  only  implosion  kinetic  energy  were 

being  converted  to  K-shell  emission.  The  comparisons  between  the 

theoretically  calculated  mass,  implosion  time,  current  at  the  time  of  peak 

kinetic  energy  generation,  kinetic-energy-per-ion,  and  K-shell  yield  and  their 

experimentally  determined  values  are  shown  in  Figures  1-5.  In  the  experiment, 

2  . 

values  of  m  were  chosen  according  to  the  relation  m  1/r  as  shown  in  Fig. 

2  ° 

1.  Note,  there  is  a  slight  deviation  from  the  mrQ  =  constant  curve  for  the 
80  yg/cm  load.  The  collapse  time  and  peak  current,  which  should  have  remained 
roughly  constant  as  rQ  and  m  were  varied,  had,  in  fact,  the  experimental 
variations  displayed  in  Figs.  2  and  3,  respectively.  Because  the  substantial 
deviations  in  the  peak  current  shown  in  Fig.  3  could  represent  a  factor  of  two 
differential  in  energy  coupling  to  the  load,  these  deviations  need  to  be  taken 
into  consideration  when  analyzing  the  data.  Fig.  4  is  a  plot  of  h  versus  2rQ, 
where  h  is  defined  as  the  kinetic-energy-per-ion  divided  by  the  minimum 
kinetic  energy-per-ion  needed  to  excite  an  ion  into  the  K-shell.  The  minimum 
kinetic-energy-per-ion  is  theoretically  predicted  to  be  about  12  keV  for  an 
aluminum  plasma.  If  this  is  a  true  prediction,  then  one  would  expect  a 
precipitous  drop  in  K-shell  emission  as  f|  goes  below  unity.  Indeed,  this  is 
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what  occurred  in  the  experiment  as  reflected  in  Fig.  5,  where  the  K-shell 
yields  as  a  function  of  array  diameters  are  graphed.  There  is  a  large 
decrease  in  K-shell  emission  which  occurs  between  the  12.5  mm  array,  which  has 
an  h  value  of  1.3,  and  the  9  mm  array,  which  has  a  value  of  q  below  unity. 
The  experiments  also  show  a  qualitative  change  in  the  K-shell  emission  region 
when  q  becomes  <  1. 

When  the  masses  imploded  in  these  experiments  are  large  enough,  the 
theory  predicts  that  the  K-shell  yield  should  be  in  an  I  scaling  regime  and 
consequently  the  yield  should  be  approximately  30%  -  36%  of  the  total  kinetic 
energy.  Therefore  any  large  variation  in  the  theoretical  yield  predictions 
would  be  related  to  variations  in  how  much  kinetic  energy  was  generated.  The 
calculation  of  kinetic  energy  is  based  on  estimates  of  the  changes  in  diode 
inductance  occurring  in  the  experiments  and  should  be  roughly  constant  for 
different  rQ  in  accordance  with  the  predictions  made  in  Fig.  5.  This  is 
interesting,  because  the  experiments  show  that  the  maximum  experimental  K- 
shell  yield,  which  occurs  for  q  near  but  >  1,  is  a  factor  of  3  larger  than 
predicted  by  kinetic  energy  conversion  theory,  i.e.  it  is  as  if  all  of  the 
imploded  kinetic  energy  had  been  converted  into  K-shell  emission  (which  is  not 
possible).  Consequently,  one  can  infer  from  the  data  that  there  is  more 
energy  being  converted  to  x-rays  than  was  originally  thought  to  occur  by 
kinetic  energy  conversion  alone.  This  is  plausible  since  the  current  remains 
on  in  the  experiment  after  the  implosion  has  occurred'  ’  ' .  Energy  is  also 
likely  to  be  trapped  in  the  diode  by  an  electrical  short  that  develops  as  the 
plasma  inductance  rapidly  increases  as  the  plasma  nears  peak  implosion  .  This 
trapped  magnetic  field  energy  is  transferred  to  the  plasma  either  ohmically  or 
inductively  and  it  ultimately  escapes  the  plasma  as  radiation.  The  time  over 
which  emission  occurs  should  be  indicative  of  the  decay  time  of  the  energy 
trapped  in  the  magnetic  field.  Because  no  further  collapse  of  the  plasma 
appears  to  be  observed  experimentally,  one  expects  that  an  anomalous  ohmic 
heating  is  occurring  and  that  this  heating  is  largest  for  q  near  1.  More  work 
is  being  done  to  analyze  these  experiments.  Note,  though,  that  the  aluminum 
calculations  did  correctly  predict  minimum  x-ray  output  from  the  pinch.  This 
is  an  extremely  encouraging  result. 


3 


We  have  been  motivated  by  these  comparisons  between  experiment  and 

theory  to  investigate  a  number  of  questions  and  speculations  regarding  their 

interpretation.  One  such  speculation  is  that  it  may  be  best  in  general  (and 

not  just  in  the  case  of  these  particular  aluminum  experiments)  to  implode  the 

maximum  amount  of  mass  that  one  can  while  still  allowing  h  to  be  in  excess  of 

unity.  There  are  three  reasons  for  this  speculation.  First,  of  course,  this 

is  what  the  comparison  of  K-shell  yield  versus  2rQ  shows  in  Fig.  5.  Secondly, 

the  theory  predicts  that  maximizing  mass  while  keeping  0  greater  than  1  will 

virtually  guarantee  that  maximum  conversion  efficiency  of  machine  energy  to  K- 

shell  emission  is  achieved.  The  third  reason  is  related  to  the  fact  that 

there  appears  to  be  experimental  evidence  to  support  the  conclusion  that  the 

experiments  never  achieve  as  large  a  density  as  predicted  by  theory.  Since  a 

large  density  as  well  as  a  long  pulse  duration  are  key  ingredients  in 

2 

determining  when  a  plasma  is  in  the  maximum  efficiency  (I  )  theoretical  K- 
shell  emission  scaling  regime,  it  makes  plausible  sense  that  if  the 
experiments  that  achieve  the  largest  density  have  the  more  massive  loads,  that 
they  are  thus  more  likely  to  produce  maximum  emission.  This  is  at  present 
only  a  reasonable  speculation.  However,  we  are  presently  in  the  process  of 
doing  a  comparison  between  experimental  and  theoretical  temperatures  and 
densities  in  order  to  quantify  any  differences  between  theory  and  experiment 
and  to  verify  the  applicability  of  the  above  speculation  to  these  particular 
experiments. 

We  have  also  begun  investigating  the  implosion  dynamics  of  nickel  wire 
arrays  using  the  same  simplifying  assumptions  as  were  used  in  the  aluminum 
calculations.  During  the  runin  phase  of  the  dynamics,  we  assumed  that  the 
current  was  linearly  rising.  In  this  case  the  final  kinetic-energy-per-ion 
that  is  generated  prior  to  assembly  on  axis  depends  only  on  the  implosion  time 
of  the  array.  Secondly,  we  assumed  that  the  current  was  shut  off  just  before 
the  plasma  collided  with  itself  on  axis.  Then,  the  dynamics  of  implosion  and 
collision  could  be  calculated  using  a  1-dimensional  magneto-hydrodynamics 
(MHD)  code  that  contains  three  essential  attributes:  (1)  it  contains  a 

structurally  complete,  state-of-the-art  atomic  model  of  nickel,  (or  aluminum 
for  the  work  described  above),  (2)  it  contains  an  extensive  nonlinear, 
nonlocal  description  of  line  photo-couplings  within  the  plasma.  (Because  of 
the  inelastic  nature  of  the  z-pinch  collision  process,  to  be  seen  shortly, 
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this  attribute  is  of  prime  importance),  and  (3)  it  consists  of  a  fully 
consistent  coupling  of  the  atomic  ionization  model  and  radiation  transport  to 
the  array  implosion  MHD. 

From  an  analysis  of  a  slug  model  description  of  the  implosion  dynamics, 

one  finds  that  two  energies  play  vital  roles  in  determining  the  radiation 

characteristics  of  the  imploded  array.  The  maximum  kinetic-energy-per-ion, 

which  is  drawn  in  Figure  6  for  the  radiating  cases  shown  in  Figures  7  and  9, 

controls  the  temperature  that  is  reached  by  the  plasma.  In  order  to  achieve 

greater  simulation  fidelity  with  plasma  radiation  sources  (PRS)  it  is 

necessary  to  radiate  from  the  K-shell  of  moderate-Z  elements  up  to  krypton. 

As  seen  in  Figure  6,  to  radiate  in  the  K-shell  of  aluminum  with  energies 

around  1.8  keV  requires  a  kinetic-energy-per-ion,  in  the  case  under  study,  of 

50  keV/ion  (4  times  the  minimum);  whereas,  to  radiate  in  the  K-shell  of 

nickel  around  8  keV  and  beyond  requires  that  the  nickel  ions  reach 

corresponding  energies  of  700  keV/ion.  If  these  ion  energies  are  reached,  the 

amount  of  radiation  actually  emitted  from  the  K-shell  depends  on  the  amount  of 

kinetic-energy-per-centimeter  generated  in  the  implosion,  i.e.  on  the  mass 

imploded.  For  aluminum  imploded  to  a  kinetic-energy-per-ion  of  ~  50  keV/ion, 

2 

the  yield  of  K-shell  x-rays  increases  at  first  as  m  at  low  masses  and  then  as 
m  at  high  masses  (greater  than  24  ygm/cm  as  seen  in  Figure  7).  This 
transition  occurs  as  a  result  of  the  collision  becoming  more  and  more 
inelastic  as  the  mass  increases  (see  Figure  8).  When  this  transition  takes 
place,  roughly  36%  of  the  kinetic-energy-per-centimeter  generated  during  the 
implosion  can  be  converted  to  K-shell  x-rays. 

As  seen  in  Figure  9,  the  results  for  aluminum  can  be  scaled  to  nickel. 
For  the  first  three  cases  shown  in  the  table,  the  nickel  reached  a  kinetic- 
energy-per-ion  of  ~  700  keV/ion  when  imploded  from  1  cm  in  33  nsec.  (If  one 
wanted  to  lengthen  this  implosion  time,  one  would  have  to  start  the  array 
implosion  from  proportionately  larger  radii,  hence  the  importance  of  studying 
screwpinch  implosions.)  As  the  mass  is  increased  from  10  to  100  to  500 
ygm/cm,  the  K-shell  yield  from  nickel  at  first  scales  as  m  and  then  begins  to 
approach  m-scaling.  The  peak  currents  needed  to  achieve  these  implosions 
scales  as  /m.  Thus,  20  MA  is  needed  to  implode  nickel  from  1  cm  in  33  ns  to 
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achieve  a  yield  of  58  kJ/cm.  (A  2  cm  length  array  would  radiate  116  kJ . )  As 
the  last  case  in  Figure  9  shows,  the  same  20  MA  can,  in  principle,  produce 
increased  yields  if  the  implosion  time  can  be  shortened. 

It  is  worth  contrasting  the  results  of  this  paper  with  the  results 
presented  in  another  paper  at  the  Simulation  Fidelity  Workshop.  In  C. 
Hartmann's  paper  on  "The  Compact  Torus  Accelerator,  A  New  X-Ray  Source  for 
High  Fidelity  Simulator",  we  learned  that  if  an  energetic  beam  of  ions, 
enclosed  in  some  electrons,  i.e.  a  plasma,  is  shot  at  a  stationary  copper  wall 
it  bounces  from  the  wall  and  generates,  in  principle,  upward  of  hundreds  of 
kilojoules  of  10  keV  x-rays.  In  nickel  z-pinch  implosions  we  showed  that,  if 
a  plasma  is  collided  and  compressed  against  itself  (compare  ring  accelerator 
concepts),  one  can  get,  in  principle,  hundreds  of  kilojoules  of  10  keV  x-rays 
if  one  controls  the  amount  of  mass  participating  in  the  emission  carefully. 
Both  concepts  require  that  the  plasma  be  focusable.  Both  have  their  problems 
in  this  regard,  and  both  have  different  options  for  flexibilitv.  For  this 
reason,  both  research  approaches  and  options  should  be  vigorously  pursued. 
However,  there  is  one  important  difference  between  the  two  approaches  that 
should  be  noted.  In  z-pinches  significantly  more  than  kinetic  energy 
conversion  has  been  observed  to  occur  in  the  experiments,  as  the  current 
continues  to  cook  the  plasma.  This  important  energy  dividend  is  not  available 
in  the  ring  accelerator. 
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PRELIMINARY  ANALYSIS  OF  ALUMINUM  2-PINCH  -  DATA  TAKEN  ON  BLACKJACK  5 


As  described  in  the  1988  DNA  Final  Report  [1]  and  in  NRL  Memo  6383  |2],  one 
can  use  theoretically  determined  line  ratios  and  yields  to  estimate  the  average  ion 
density  and  electron  temperature  reached  during  peak  compression  in  a  particular 
experiment.  In  previous  work,  ve  looked  at  titanium  vire  experiments  performed  on 
Blackjack  5.  Ue  have  now  taken  a  preliminary  look  at  several  aluminum  vire 
experiments  recently  done  on  Blackjack  5. 

As  with  the  theoretical  titanium  database  used  previously,  the  aluminum 
database  is  composed  of  a  series  of  collisional  radiative  equilibrium  (CRE) 
calculations  done  for  a  steady-state,  homogeneous  cylindrical  plasma  at  different 
ion  densities,  electron  temperatures  and  radii.  From  these  CRE  calculations  ve  have 
derived  the  dependence  of  the  radiated  power  per  centimeter,  the  He-a  to  Ly-a  line 
ratio,  the  He-(3  to  Ly-a  line  ratio,  and  the  He-3  to  Ly-g  line  ratio  on  the  ion 
density  and  electron  temperature.  Given  the  size  of  the  emitting  region,  these  four 
values  can  then  be  compared  to  the  experimentally  determined  average  K-shell  energy 
production  rate  and  the  relative  line  strengths. 

2.  Experimental  Data 

Seven  of  the  filtered  SXR  series  aluminum  vire  array  shots,  recently  performed 
by  Nick  Loter  on  BLACKJACK  5  at  Maxwell  Laboratories,  have  been  analyzed.  These 
shots  were  all  filtered  by  the  addition  of  1.5  mils  of  kapton  and  all  except  one 
shot,  shot  2137,  contained  some  magnesium.  The  mass  loading  ranged  from  about  80  to 
400  ug/cm,  while  the  resulting  K-shell  yields  ranged  from  around  7  to  40  KJ.  Table 
1  summarizes  the  experimental  parameters. 

TABLE  1.  BLACKJACK  5  source  performance  during  aluminum  vire  SXR  series. 


Shot 

Mass 

K-shell 

K-shell 

Emitting 

Plasma 

No. 

Loading 

FVHM 

Yield 

Radius 

Length 

(ug/cm) 

(ns) 

(KJ) 

(cm) 

(cm) 

2111 

53 

30 

8.6 

0.065 

3.1 

2122 

118 

37 

20.6 

0.045 

3.1 

2125 

161 

57 

32.3 

0.050 

3.1 

2126 

210 

49 

20.7 

0.035 

3.1 

2127 

210 

51 

32.9 

0.050 

3.1 

2132 

210 

48 

29.6 

0.050 

3. 5(2.1) 

2137 

397 

60 

15.9 

0.100 

1.7 

Based  on  a  table  and  other  information  provided  by  Nick  Loter,  Maxwell  Labs. 
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Both  the  emitting  radius  and  plasma  length  were  estimated  from  pinhole 
pictures.  Due  to  the  resolution  of  the  pinhole  pictures,  accurate  measurements 
were  difficult,  therefore,  the  numbers  quoted  in  Table  1  are  approximate.  The 
actual  emitting  radius  can  vary  considerably  along  the  length  of  the  plasma 
column.  Indeed  in  some  cases,  the  emission  region  can  pinch  off  entirely 
forming  almost  spherical  blobs. 

From  table  1  and  the  spectra,  the  average  K-shell  energy  production  rate 
and  the  line  ratios  can  be  determined.  These  values  are  shown  in  Table  2. 
Note  that  the  K-shell  yields  in  Table  1  include  a  small  amount  of  energy  from 
the  the  magnesium  lines.  To  account  for  this  we  have  decreased  the  K-shell 
yields  for  all  shots  except  shot  2127  by  5%.  As  with  the  emitting  radius  and 
length,  there  is  some  as  yet  unknown  amount  of  uncertainty  in  these  values. 

TABLE  2.  K-shell  Production  Rate  and  Line  Ratios  for  aluminum  shots. 


Shot 

No 

J/ns/cm 

He-a/Ly-a 

He-8/Ly-a 

He-B/Ly-B 

2111 

88 

0.66 

0.12 

1.11 

2122 

171 

0.68 

0.11 

0.97 

2125 

174 

0.86 

0.13 

1.41 

2126 

129 

1.47 

0.25 

0.57 

2127 

198 

1.26 

0.11 

3.16 

2132 

167(288) 

0.68 

0.14 

2.10 

2137 

156 

1.38 

0.08 

2.06 

Analysis  of  Experimental  Data 

A  fixed  value  of  radiated  power  per  centimeter  will  map  out  a  single 
contour  in  (N^,Te)  space,  where  N^  is  the  average  ion  density  and  Tg  the 
average  electron  temperature.  Similarly,  the  three  line  ratios  will  map  out 
three  different  contours.  The  intersection  of  the  contour  for  a  fixed  radiated 
power  per  centimeter  with  any  of  the  contours  for  a  fixed  line  ratio  will  give 
a  single  ion  density  and  electron  temperature.  One  can,  therefore,  obtain  the 
average  conditions  needed  to  produce  each  of  the  line  ratios  seen  in  the 
experimental  spectrum.  By  averaging  all  three  temperatures  or  densities 
together,  one  would  obtain  an  overall  average  temperature  and  density  necessary 
to  produce  the  K-shell  emission  seen.  Figure  1  illustrates  this  for  filtered 
shot  2125. 
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A  separate  estimate  of  the  electron  temperature  can  be  made  from  the 
continuum.  All  four  temperature  estimates  are  shown  in  Table  3,  while  the 
three  ion  densities  are  listed  in  Table  4.  The  values  quoted  depend  on  the 
line  ratios,  size  of  the  emitting  region  and  yield.  Any  uncertainties  in  these 
quantities  will  result  uncertainty  in  the  estimate  of  the  average  temperature 
and  ion  densities  reached  during  peak  compression. 

TABLE  3.  Estimates  of  the  Electron  Temperature  (in  eV). 

Shot  From  From  Intersection  of  J/ns/cm  with: 

No  Continuum  He-a/Ly-a  He-6/Ly-a  He-6/Ly-g 


2111 

770 

430 

2122 

800 

570 

500 

2125 

520 

490 

460 

2126 

435 

380 

410 

2127 

400 

490 

2132 

640 

570(540) 

450(480) 

2137 

520 

450 

TABLE  4.  Estimates  of  the  Ion  Density  (cm  )  from  Intersection  of  J/ns/cm  with 
the  three  line  ratios. 


with 

He-a/Ly-a 

2.9x10^ 

1. 4x102- 

1.5x10^ 
2.6x10,^ 
2.4xl(TU  , 
1.2(2.21x10' 


2137  5.0x10 

Conclusions 


He-6/Ly-a 

6.2x10:^ 

1. 7x10^0 

i.7xio;^ 

2.4x10^ 

1.7x10^  , 

1.7(2.51x10' 

3.9xKT° 


He-g/Ly-3 

6.2x10^ 

2.1x10^ 

2.2x10^ 

1.7xl0ZU 

3.4(5.31x10' 

9.0x10 


Using  the  ion  density  obtained  from  the  intersection  of  the  radiated 
power  per  centimeter  (J/ns/cm)  with  the  He-a  to  Ly-a  line  ratio,  one  can 
estimate  the  amount  of  mass  that  is  radiating.  This  quantity  is  shown  in  Table 
5.  The  last  column  gives  the  fraction  of  the  initial  mass  that  radiated. 

An  estimate  of  how  sensitive  the  average  electron  temperature  and  ion 
density  are  to  uncertainty  in  the  size  of  the  emitting  region,  the  yield  and 
the  line  ratios  still  needs  to  be  done. 
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TABLE  5.  Emitting  Mass  versus  Initial  Mass. 


Shot 

Initial 

Emi t  ting 

Fraction 

No 

Mass 

Mass 

Emi 1 1 ing 

(Ug/cm) 

(Ug/cm) 

% 

2111 

53 

17 

33 

2122 

118 

40 

34 

2125 

161 

53 

33 

2126 

210 

45 

21 

2127 

210 

84 

40 

2132 

210 

42(77) 

20(37) 

2137 

397 

70 

18 
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n.  Atomic  Data  Aquisition 


Complete  data  bases  of  collision  strengths  for  the  hydrogen-like  to  neon-like  ionization 
stages  of  iron  and  selenium  have  been  finished.  Our  next  goal  in  the  area  of  data  bases  is  to 
complete  the  ionization  stages  of  H-like  to  Ne-like  for  argon.  The  data  base  will  include 
energy  levels,  transition  probablities,  photoionization  cross  sections,  ionization  and 
excitation  collision  strengths  and  rates.  This  data  will  allow  us  to  use  the  Z-scaling 
software  that  is  under  development  to  determine  oscillator  and  collision  strengths  for  18  < 
Z  <  36.  Otherwise,  it  would  be  prohibitive  in  time  and  storage  either  to  calculate  or  store 
all  the  necessary  oscillator  and  collision  strength  data  for  this  large  number  of  elements  and 
ionization  stages.  It  is  sufficient  to  provide  fits  to  a  much  smaller  data  base  that  is  accurate 
enough  for  modeling  and  transport  calculations  moreover  there  is  a  built  in  reliability  to  this 
procudure  for  producing  data. 

The  software  used  to  discover  the  systemic  trends  which  exist  for  atomic  oscillator 
strengths  has  been  developed  and  implemented  successfully.  The  electric  dipole  oscillator 
strengths,  computed  using  R.  D.  Cowan’s  atomic  structure  package,  show  smooth 
variations  of  the  calculated  values  as  a  function  of  the  nuclear  charge,  Z.  The  oscillator 
strength  fit  is  of  the  form 


f,z' 


+  f2z-2 


where  fQ  is  the  hydrogenic  oscillator  strength.  The  helium-like  and  neon-like  models  have 
been  completed. 

The  results  from  Cowan's  Program  RCG  Mod  9  fit  remarkably  well  into  the  systematic 
trends  for  both  of  the  completed  models.  Table  1  gives  the  details  of  each  model.  Figure 
1  is  a  plot  of  f-number  versus  nuclear  charge  for  the  He-like  model.  Transitions  A  and  L 
show  smooth  variations  with  nuclear  charge.  As  was  expected,  these  two  transitions  are 
simply  a  multiplicative  factor  of  the  hydrogenic  values.  Curve  K,  as  well  as  all  other  An=0 
transitions,  are  distincitive  in  curvature  and  exhibit  more  pronounced  1/Z  behaviour.  This 
is  a  result  of  the  hydrogenic  value,  fQ  ,  being  quite  small.  In  hydrogen,  An=0  transitions 

are  degenerate  which  means  the  oscillator  strengths  are  zero.  Figure  2  shows  two 
transitions  from  the  Ne-like  model. 

The  oscillator  strength  models  that  we  are  building  have  the  advantage  of  being 
completely  self-consistent.  In  other  words,  the  random  errors  which  have  plagued  similar 
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calculations  and  are  a  result  of  gathering  avialable  data  from  various  sources  have  been 
completely  eliminated. 


Table  1 :  Oscillator  Strength  Models 


Isoelectronic 

Number  of 

Number  of 

Transition 

Data 

Sequence 

Energy 

Transition 

Arrays 

Base 

Modeled 

Levels 

Probabilities 

He 

38 

80 

Is2  -ls1nl 

Ne  IX, 

ls!nl  -  lsVl' 

Ar  XVII, 

Fe  XXV, 

Se  XXXIII 

Ne 

51 

83 

2s22p6  -  2s22p5nl 

Ar  IX, 

2s22p6  -  2s12p6nl 

Ti  XIII, 

2s22p5nl  -  2sz2pVl' 

Fe  XVII, 

2s22p5nl  -  2s12p6nl 
2s12p6nl  -  2s12p6nT 

SeXXV 

Finally,  the  development  of  software  to  fit  collision  strengths  along  isoelectronic 
sequences  is  well  underway.  The  form  that  we  are  attempting  to  fit  for  the  collision 
strength  data  is  a  two  parameter  fit  in  energy  and  nuclear  charge. 

Q(X)  =  F1(Z)(A  +  BX_1+  CX'2)  +  F2(Z)  ln(X) 

2 

where  Fn(Z)  =  (Z  -  Bn)  are  screening  factors  and  X  is  energy  in  threshold  units. 

The  first  stage  in  the  process  is  an  energy  fit  for  each  of  the  transitions.  A  simple 
procedure  for  the  energy  fitting  to  the  data  has  been  implemented  for  the  Ne-like 
isoelectronic  sequences  of  iron  and  selenium.  A  comparison  of  the  two  data  sets  provided 
a  check  for  the  relative  correctness  of  the  data.  Because  the  matrix  element  is  proportional 
to  Z  ,  the  collision  strength  magnitude  should  get  smaller  as  Z  increases.  This  fact  led  us 
to  discover  a  serious  error  for  a  few  of  the  transitions  in  the  iron  data  base.  This  is 
illustrated  in  figure  3.  The  curve  labeled  C  is  the  old  iron  data  and  is  much  lower  than  the 
selenium  data  for  this  transition  (curve  B).  Notice  after  correcting  the  mistake,  the  new 
iron  data  (A)  and  the  selenium  data  are  in  the  correct  relative  order.  With  the  Ne-like  iron 
corrected,  we  can  proceed  to  the  next  stage  in  the  scaling  process. 
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DIELECTRONIC  RECOMBINATION  RATE  COEFFICIENTS  FOR  F-LIKE  IONS  AND  Z-SCALING 


OF  AUTOIONIZATION  AND  DIELECTRONIC  RECOMBINATION  RATES 

Dielectronic  recombination  (DR)  is  an  important  recombination  process  for 
plasma  diagnostics  of  atomic  processes  in  high  temperature  and  moderate 
density  plasmas.  The  dielectronic  satellite  lines  observed  in  the  emission 
spectra  as  a  result  of  these  recombinations  can  be  used  for  such  plasma 
diagnostics.  The  models  used  to  describe  high  temperature  plasmas  consist 
of  rate  coefficients  for  atomic  processes  in  the  ionization  balance  or 
transport  equations.  Semi-empirical  formulas  predict  the  total  effective 
ground-to-ground  DR  rate  coefficients  at  best  within  a  factor  of  two.  On 
the  other  hand,  the  calculations  being  carried  out  here  for  the  first  time 
will  provide  accurate  DR  rates  into  specific  excited  states  (i.e.  into  the 
states  that  radiate)  and  from  which  the  highly  approximate  total  effective 
rate  can  be  obtained,  if  needed,  after  summing  the  state  specific  rates. 
More  improved  calculations  of  DR  and  accurate  calculations  of  other  rate 
coefficients  are  thus  extremely  important  for  their  utility  as  proper 
plasma  diagnostics  and  for  accurate  calculations  of  X-ray  yields  in  the  L- 
shell . 

Dielectronic  recombination  is  a  two  step 
process.  In  the  first  step  ,  there  is  an  initial  capture  of  an  electron  by 
an  ion  X+z  into  a  doubly  excited  state  (X+^Z-^)*  which  stabilizes  by 
subsequent  radiative  decay  to  singly  excited  states  (X+^z  ^)  or  by 
autoionization  to  the  initial  state  X+z  or  any  excited  stated  (X+z)  of  the 
initial  ion.  In  order  to  calculate  the  total  DR  rate  from  one  ion  stage  to 
the  next  or  even  the  total  DR  rates  to  a  specific  singly  excited  state  of 
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the  next  ion  stage,  one  has  to  calculate  from  a  large  number  of  doubly 
excited  state  all  possible  radiative  decay  rates  as  well  as  all 
energetically  allowed  autoionization  rates  .  This  amounts  to  calculating  an 
enormous  amount  of  radiative  and  autoionization  rate  data.  However,  for  an 
accurate  calculation  of  the  DR  rates,  a  detailed  level  accounting  scheme, 
with  extrapolations  or  configuration  averaging  only  when  sufficient 
convergence  is  present,  is  crucial. 

In  this  years  work  we  have  completed  the 
calculation  of  DR  for  ions  of  the  flourine  isoelectronic  sequence  covering 
ions  with  Z=  18,22,26  and  34.  Ve  have  also  calculated  DR  rates  from  neon 
like  ions  and  we  are  thus  in  the  process  of  developing  scaling  laws  for 
these  rates  throughout  the  K-  and  L-shell  ionization  stages.  We  calculate 
DR  rates  by  calculating  autoionization  rates  and  radiative  rates  from 
intermediate  resonance  states  using  a  Hartree-Fock  calculation  with 
relativistic  correction  (  HFR  method  of  R.  D.  Cowan  Approximate 
relativistic  corrections  in  the  radial  wave  functions  as  well  as  in  the 
energies  are  included  in  this  calculation.  Thus  calculation  of  DR  rates 
involve  calculations  and  organizations  of  a  range  of  atomic  quantities 
such  as  energies  of  the  doubly  excited  and  singly  excited  states, 
autoionization  and  radiative  rates  from  each  doubly  excited  states  to  each 
recombining  states  and  final  stabilizing  states  respectively.  We  thus 
developed  computer  codes  to  extract  relevant  quantities  from  the  output 
files  of  this  structure  code  calculations  and  stored  and  organized  them 
into  different  data  files  so  as  that  we  can  use  them  not  just  to  obtain  DR 
rates,  but  even  configuration  averaged  autoionization  rates  ,  radiative 
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decay  rates,  and  resonance  excitation  rates  from  a  particular  state  to 
another. 

This  amounted  to  a  large  amount  of  bookkeeping  and  organizational  work 
(which  had  to  be  done  only  once  )  since  depending  on  the  complexity  of  the 
ion  there  could  be  a  very  large  number  of  doubly  and  singly  excited  states. 
In  calculating  the  DR  rates  we  are  thus  obtaining  data  which  will  give  us 
other  important  rates  besides  DR  rates  such  as  resonance  excitation  rates 
as  mentioned  before  or  quantites  needed  for  dielctronic  satellite  spectra 
calculations. 


Ue  expect  that  this  method  of  calculation  will  give  more  accurate  results 

than  previously  calculated  DR  rates  using  single  conf iguaration  non- 

relativistic  wave  functions  or  using  a  simple  angular  momentum  averaged 

procedure  by  averaging  over  many  intermediate  states.  On  the  other  hand  it 

DR 

is  less  time  consuming  and  less  complex  than  calculating  a  using  fully 

relativistic  multiconfiguration  Dirac-Fock  (MCDF)  method.  Thus  for 

calculating  autoionization  or  DR  rates  for  moderate  Z  ions,  our  method  of 
DR 

calculating  a  is  a  very  good  compromise  between  accuracy  and  computer 
time. 

DR 

The  DR  rate  coefficient  a  from  an  initial 
state  |i>  into  a  specific  stabilizing  bound  state  |k>  through  a  particular 
autoionization  state  |j>  is  given  by  : 


D*  .  . 

cK  (  i,  l<  k} 


a3; 


o; 
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for  a  Maxwellian  distribution  of  electrons  ;  where  ang  g^  are  the 

statistical  weights  of  the  initial  and  doubly  excited  states  respectively, 

T  is  the  electron  temperature  and  is  the  energy  of  the  free  recombining 

s  r 

electron  in  Rydbergs.  A  ^  and  A  ^  are  the  autoionization  and  radiative 
rates  from  state  |j>  to  states  |i>  and  |k>  respectively. 

The  Auger  transition  probability  for  each  doubly 
excited  state  is  calculated  from  the  perturbation  expression 


A..a  =(  2n/h  )  |<j|l/r12|i>|2 


where  |j>  and  |i>  are  intermediate-coupling  wave  functions  and  the 
radiative  electric  dipole  transition  probability  from  state  |j>  to  |k>  is 
given  by 


Ajkr  =  (  2it/2Jj  +  l  )  |<k|T|j>|2  (3) 

where  <k|T|j>  is  the  reduced  matrix  element.  The  transition  matrix  elements 

in  Eqs.  (2)  and  (3)  can  be  expressed  as  sums  of  terms  involving  the  product 

of  angular  integrals  and  radial  integrals.  These  matrix  elements  were 

calculated  using  Cowan's  atomic  stucture  code^.(RCG  Mod  9) 

(2  3) 

Most  published  works  ’  report  only  the  total 

pU  . 

dielectronic  recombination  rates  T where  the  summation  over 

<«*,  j,  K 

m  extends  over  all  levels  of  the  ground  configuration  |i>  of  the  ion,  as  a 
function  of  temperature  and  the  sum  extends  over  all  the  singly  excited 
bound  final  states  as  well.  However,  separate  bookkeeping  and  presentation 
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from  state 


of  partial  atomic  data  such  as  the  autoionization  rates 
|j>  to  state  |i>,  radiative  rates  A^1"  from  state  |j>  to  all  singly 
excited  levels  |k>  below  the  ionization  limit  and  dielectronic 

recombination  rates  initial  state  |i>  to  specific  final  state 

.  i 

|k>  are  extremely  important  for  various  reasons.  Vith  only  the  total  rate 
coefficients,  it  is  not  possible  to  study  density  effects  on  these  rates 
and  to  calculate  DR  rates  accurately  at  densities  of  PRS  (Plasma  Radiation 
Source)  interest,  dielectronic  recombination  rates  to  these  specific  singly 
excited  bound  states  are  needed.  Therefore,  considerable  effort  was 

expended  t'1  set  up  the  computer  bookkeeping  to  isolate  these  rates  and  to 
define  and  calculate  the  necessary  data  files  to  make  this  data  available 
to  other  of  the  ionization  calculations  that  are  to  be  carried  out  in  this 


DNA  program. 


To  calculate  Dielectronic  recombination  from 


F-like  ion  to  a  Ne-like  ion  the  initial  continum  states  before 
recombination  are  taken  to  be: 

ls22s22p5  (2P)  and  ls22s2p6  (2S) 

The  autoionization  states  for  An*0  , where  n  refers  to  principal  quantum 
number  of  the  excited  electron  are  : 

2s22p431n'l'  ,  2s2p531n'l'  and  2s°2p631n'l' 

These  states  stabilize  by  dipole  allowed  radiative  transitions  to  the 


following  states: 


2s22p^n'l'  and  2s2p^n'l' 


where  n'  and  1'  represent  the  quantum  numbers  of  the  captured  electron. 
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For  An=0  transitions,  the  autoionization  states  are  of  the  form  2s2p^n'l' 

2  5  6 

and  these  can  stabilize  by  radiative  decay  to  2s  2p  n"l",  2s2p  n''l''  and 

2  (j 

to  the  ground  state  2s  2p  of  Ne-like  ion.  These  transitions  become 

energetically  possible  for  high  Rydberg  states,  specifically  for  higher  Z 

ions  such  as  selenium.  Additional  autoionization  to  excited  states  of  the 

2  4  5 

recombining  ion  included  2s  2p  31  and  2s2p  31  states  whenever  they  were 
energetically  possible. 

In  this  present  calculation,  we  neglected  the  cascade 
effect,  i.e,  radiative  transitions  leading  to  other  autoionization  states 
since  it  is  assumed  that  its  contribution  in  reducing  the  DR  coefficients 
is  very  small. 

r*  a 

For  An*0,  we  calculated  the  radiation  rates  A  and  autoionization  rates  A  , 
DR 

i.e,  a  explicitly  for  n<  10  and  l'<  6.  For  states  with  n>10  ,  the  DR 

3 

coefficients  were  obtained  by  using  the  1/n  extrapolation  to  these 
coefficients.  For  An=0,  as  mentioned  before,  autoionization  transitions 
become  energetically  possible  for  n>  nQ  where  nQ  varies  from  ion  to  ion. 
For  each  ion  the  value  of  nQ  was  determined  by  computing  the  energy  levels 
of  the  autoionizing  state.  We  calculated  the  transition  rates  for  n  <  n<  15 

q-  “ 

and  1<  8  for  these  An=0  transitions.  For  states  with  n  >15,  the  radiative 
rates  remained  constant  where  as  the  autoionization  rates  were  extrapolated 

3 

by  the  usual  1/n  extrapolation  procedure.  For  An=0,  radiative  decay  from 
the  Rydberg  electron  is  found  to  contribute  significantly  to  the  DR  rate 
coefficients.  For  the  autoionization  states  of  the  2p^31nl'  configurations, 
these  states  can  autoionize  to  the  excited  states  of  the  initial  ion  as 
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mentioned  before.  Ue  have  also  included  these  autoionization  channels  in 


our  calculations. 


These  rates  were  calculated  using  the 
structure  code  of  R.D.Cowan^.  The  calculations  were  done  in  the  Hartree- 
Fock  approximation.  Relativistic  terms  were  included  in  the  differential 
equation  to  give  relativistic  correction  to  the  radial  wave  functions  as 
well  to  obtain  improved  energies.  The  continum  wave  functions  were 
calculated  in  the  distorted  wave  approximation.  All  calculations  were 
carried  out  in  intermediate  coupling  and  single  configuration 
approximation.  Thus  configuration  interaction  between  autoionizing  states 
were  neglected  in  this  calculation. 

DR  rate  coefficients  for  the  F-like  ions  were 
calculated  for  ions  with  Z=18,26  and  34.  Ue  have  also  completed  a  similar 
calculation  for  F-like  titanium  with  Z=22.  However,  the  results  for  Ti  is 
not  included  in  this  report.  These  rates  were  calculated  and  new 
bookkeeping  procedures  were  set  up  using  the  atomic  structure  code  of  R.  D. 
Cowan^^.  Data  files  were  created  for  each  ion  with  a  specific  structure 
for  calculating  and  storing  the  data.  Ue  have  calculated  and  stored  the 
energies  of  the  doubly  and  singly  excited  Ne-like  states  with  respect  to 
the  ground  state  of  F-like  state.  The  energy  level  structure  of  the  Ne-like 
autoionizing  and  singly  excited  states  are  built  such  that  high  Rydberg 
states  with  closely  lying  energies  are  grouped  together  and  identified  as 
one  level.  In  this  present  caculation,  we  include  all  autoionization  states 
with  Rydberg  electron  having  principal  quantum  number  upto  n=10  or  15  for 
An *0  and  An=0  respectively.  Ue  included  all  possible  Auger  decay  and 
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radiative  decay  channels  including  autoionization  to  excited  states. 
Radiative  decay  to  other  autoionization  states,  i.e,  cascade  effects  were 
not  included  in  this  calculation.  It  has  been  seen  from  other  DR 
calculations  that  the  contribution  of  this  cascade  effecft  is  only  about  3- 
4£  of  the  total  DR  rate.  Ve  calculated  and  stored  from  each  doubly  excited 
state  ,  all  possible  autoionization  rates,  and  all  possible  radiative  and 
dielectronic  recombination  rate  to  each  specific  singly  excited  state. 
After  calculating  the  individual  autoionization  rates  and  radiative  rates 
ve  calclate  and  store  the  branching  ratio 

>  <>r>  K 

2  fa  -  2  fa 

relevant  for  calculating  the  DR  rate  coefficient.  The  DR  rate  coefficient 
is  then  calculated  from  each  doubly  excited  state  to  each  singly  excited 
state  for  each  electron  temperature  in  the  range  0.01<T<10  keV  by  using  Eq. 
(1).  The  partial  DR  rates  are  then  summed  over  the  final  singly  excited 
states  to  obtain  the  total  DR  rate  from  F-like  ground  state  to  Ne-like 
state. 

In  Figs.  1-3  ve  present  these  DR  rate  coefficients  for 
Ar^+,  Fe^+,  and  Se^+  for  the  (An*0),  ^P  (An=0)  and  initial  states 
as  a  function  of  electron  temperature.  It  can  be  seen  from  the  graphs  that 
the  peak  values  of  the  rate  coefficients  is  at  higher  temperatures  for  ions 
vith  higher  Z  values.  From  Fig.  2  ve  also  notice  that  for  An=0  transitions, 
the  DR  rates  are  very  large  at  very  lov  temperature  and  becomes 
insignificant  at  higher  temperatures.  This  is  because  the  doubly  excited 


A  o’.j.k)  =  1 3; 
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states  for  these  transitions  are  lying  just  above  the  initial  ground 
state. 

For  An*0  ,  the  low  lying  doubly  excited  states 
i.e, states  with  Rydberg  electron  having  low  principal  quantum  number 
(31nl',n<4)  contribute  30-60%  to  the  total  rates  for  Ar^+  through  Se^+ 
respectively.  Thus  for  higher  Z  ions,  the  contribution  of  low-lying  states 
is  much  more  significant  compared  to  those  with  high  Rydberg  states.  The 
partial  DR  rate  coefficients  from  31nl'  configurations  are  compared  in  Fig. 
4  for  Z=18  through  34.  It  is  seen  from  this  figure  that  the  Mfc-h  n  values 
become  important  for  higher  temperature. 

The  dielectronic  recombination  for  the  F-like  ions  have  been 

investigated  by  different  authors  for  several  ions  .  Roszmanv  '  has  done  a 

non-relativistic  Hartree-Fock  calculations  where  the  orbital  energies, 

autoionization  and  radiative  rates  are  calculated  in  the  single 

(3) 

configuration,  LS-coupled,  distorted  wave  approximation.  Chen  on  the 
other  hand  has  calculated  DR  rate  for  F-like  ions  using  a  fully 
relativistic  multiconfiguration  Dirac-Fock  model(MCDF).  They  have  taken 
into  account  configuration  interaction  among  states  in  the  same  complex. 
Figs.  5-7  present  comparisons  of  our  results  with  the  other  two 
calculations.  In  Fig.  5  ,  we  compare  our  results  for  Ar^+  and  Fe^+  (^P, 
An*0)  with  the  nonrelativistic  HF  calculation  of  Ref.  2  .  We  see  that  our 
rates  are  always  higher  (  by  35%  at  the  maximum  value  at  0.2  and  0.4  keV 
for  Ar^+  and  Fe^+  respectively).  In  Fig.  6  we  compare  our  results  for 
Fe^+  and  Se^"*+  (^P,  An#0)  with  MCDF  calculation  of  Chen^^.  Our  results 
are  very  close  to  the  results  of  Ref.  3.  (higher  by  less  than  8%  at  the 
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1 7  +  25+ 

maximum  value  for  both  Fe1  and  Se'1  ).  All  results  of  ref.  are  the 

angular  momentum  average  results.  Even  though  ve  expected  to  be  closer  to 

the  nonrelativistic  HF  results  of  Roszmanv  compared  to  the  relativistic 

results  of  Chen^  ,  the  much  closer  agreement  with  Chen^  suggest  that 

relativity  may  not  play  a  significant  part  in  the  DR  rate  coefficients.  The 

differecne  between  our  results  and  those  of  Roszman's  could  be  due  to 

the  sensitivity  of  the  differences  in  the  atomic  models  used  in  each  case. 

2 

In  Fig.  7  we  compare  the  total  DR  rate  coefficients  for  the  initial  P 
state  of  Fe^  +  .  Again,  results  of  this  calculation  is  closer  to  those  of 
Chen^^  than  the  non-relativistic  calculation  of  Roszman^) 
In  Figs.  8-12  we  present  graphs  of  energies, autoionization  rates, 
radiative  rates  and  DR  rate  coefficients  for  several  doubly  excited  states 
as  a  function  of  Z.  Thus  using  these  graphs,  energies  of  doubly  excited 
states  ,  autoionization  rates,  radiative  rates  and  DR  rate  coefficients  can 

l 

be  found  for  the  particular  Z  of  interest.  In  fig.  8  we  present  the 
energies  of  the  3131'  doubly  excited  states  as  a  function  of  the  atomic 
number  Z.  In  Fig.  9  we  show  the  Z-scaling  of  the  oscillator  strength  for 
the  3131'  configurations  for  n=3-»2  radiative  decay  rates.  The  oscillator 
strength  can  be  written  as 

t  =  <  gu  /  gl  ^  <mc/8rt  e  )  *  Aul 

where  gy  and  g^  are  the  statistical  weights  of  the  upper  and  lower  states 

ra/J 

respectively,  A  is  the  wave-length  of  transition  and  Au^  is  the 
radiative  decay  rate.  The  oscillator  strength  scales  like 
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,rad  ,n  ,n2 

f  =  a0  +  al/Z  +  a2/Z 

The  oscillator  strength  is  mostly  independent  of  Z  whereas  the  radiative 
4 

rate  scales  as  Z  .  It  can  be  also  justified  that  the  autoionization  rates 
are  independent  of  Z.  Ve  therefore  define  a  f-number  for  the  autoionization 
rate  equivalent  to  the  oscillator  strength  for  radiative  decay. 

fai  =  (  gu  /  g1  )  (mc/8n2e2)  aQ2  Aulai 

=  bQ  +  b:/Z  +  b2/Z2 
3  i 

where  is  the  autoionization  transition  rate  for  transition  from  state 

u  to  state  1.  In  Fig.  10  and  Fig.  11  we  show  the  Z-dependence  of ^^ese 

autoionization  f -numbers  for  autoionization  from  3131'  doubly  excited 

2  2 

states  to  the  first  excited  (  S)  and  ground  (  P)  states  respectively.  Thus 

using  Cowan's  code  ,  we  have  found  the  scaling  relations  with  atomic  number 

of  energies,  autoionization  rates,  and  radiative  rates  from  doubly  excited 

states  in  Ne-like  ionization  stage.  From  these  scaling  relations,  those  for 

the  Dielectronic  recombination  rates  into  the  neon  like  singly  excited 

states  can  be  found  for  any  temperature.  In  Fig.  12  we  present  the  Z- 

2  5 

dependence  of  the  DR  rates  from  ground  P  state  to  2p  3s  excited  state  for 
temperatures  T=  100,  200,  300,  and  600  eV.  In  developing  the  Z-scaling  of 
autoionization,  radiative  and  DR  rates,  only  3131'  autoionization  doubly 
excited  and  n=3  singly  excited  Ne-like  states  have  been  considered.  It  is 
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expected  that  the  higher  autoionizing  states  (  31nl'  with  n>  4  and  nlnl' 
with  n,n'>  4)  and  singly  excited  states  with  n>  4  will  also  scale 
accordingly. 

Because  of  strong  heating  of  a  plasma,  the  distribution 

of  atoms  among  the  various  ionization  stages  is  no  longer  best  represented 

by  the  Maxwellian  electron  distribution  function.  When  the  excitation  and 

ionization  rate  coefficients  are  integrated  over  non-Maxwellian  rather  than 

the  usual  Maxwellian  electron  distribution  function,  they  do  not  detail 

balance  with  the  deexcitation  and  recombination  (  including  dielectronic 

recobmination  )  rate  coefficients.  This  leads  to  significant  shifts  in 

ionization  abundances  as  a  function  of  temperature  and  in  the  excited  state 

DR 

population  distribution.  Thus  using  the  branching  ratio  A  obtained  for 
this  F-like  calculation  and  integrating  it  over  Maxwellian  as  well  as  non- 
Maxwellian  distribution  function,  we  found  important  effects  of  these  rates 
on  the  excited  state  population  of  F-like  ions. 

We  will  soon  complete  our  Z-scaling  calculation  by 
scaling  the  energies  and  autoionization  and  radiative  transition  rates  for 
doubly  excited  states  with  high  Rydberg  electrons.  With  the  use  of  such 
scaling  laws,  it  will  be  possible  to  carry  out  more  accurate  analyses  of 
the  feasibility  of  different  elements  in  the  isoelectronic  sequence  as 
plasma  sources  for  Z-pinches.  Equally  important,  scaling  will  allow  us  to 
obtain  these  rates  throughout  the  portion  of  the  periodic  table  of  interest 
to  DNA  Z-pinches  with  a  great  savjng  in  the  amount  of  theoretical  effort 
expended  than  would  otherwise  be  needed  if  they  were  calculated  directly  as 
was  done  with  Ar,  Ti,  Fe,  and  Se. 
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DIELECTRONIC  RECOMBINATION  RATE  COEFFICIENT  (1Q-12c 


T(keV) 
Figure  1 
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DiELECTRONIC  RECOMBINATION  RATE  COEFFICIENT  (10'12cm3/S) 


T(keV) 
Figure  2 
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DIELECTRONIC  RECOMBINATION  RATE  COEFFICIENT  (1Q-12cm3/S) 


T(keV) 
Figure  3 
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DIELECTRONIC  RECOMBINATION  RATE  COEFFICIENT  (10-12cm3/S) 
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Figure  4 
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DIELECTRONIC  RECOMBINATION  RATE  COEFFICIENT  (10-12cm3/S) 
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DIELECTRONIC  RECOMBINATION  RATE  COEFFICIENT  (10" 
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DIELECTRONIC  RECOMBINATION  RATE  COEFFICIENT  (1Cr12cm3/S) 


T(keV) 


Figure  7 
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doubly  excllod  state  energies 


Z 

Figure  8 
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AUTOIONIZATION  TO  FIRST  EXCITED  (ZS)  STATE 


AUTOIONIZATION  TO  GROUND  OP)  STATE 
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III.  RADIATIVE  COLLAPSE 


One  of  the  important  objectives  of  the  Radiation  Hydrodynamics 
Branch  is  to  promote  understanding,  interpretation,  and  prediction 
of  the  dynamics  and  radiation  from  a  Z-pinch  implosion  based  upon 
numerical  simulations.  Such  studies  can  lead  to  new  methods  and 
approaches  for  improved  coupling  betveen  the  pulse  power  machine  and 
the  pinch  load,  as  well  as  enhancement  in  the  radiative  yields  from 
the  implosion.  The  ideal  configuration  of  the  machine  and  load 
would  lead  to  a  radiative  collapse  wherein  a  large  fraction  of  the 
energy  delivered  to  the  load  is  returned  as  radiation  in  the  keV 
range,  while  only  a  small  fraction  is  lost  to  turbulent  kinetic 
motions.  For  example,  a  3  MAmp  machine  could  ideally  yield  over  a 
hundred  kJoules  in  keV  radiation  from  the  radiative  collapse  of  an 
aluminum  implosion.  The  central  question  is  how  such  an  ideal 
machine-load  configuration  which  exists  in  principle  can  be  realized 
in  practice.  To  this  end  we  have  been  studying  the  Z-pinch 
implosion  of  wire  arrays  coupled  to  a  realistic  circuit  which  models 
the  pulse  power  machine. 

Unfortunately,  it  is  known  that  ab  initio  numerical  simulations 
of  PRS  Z-pinch  loads  with  radiation  hydrodynamic  codes  do  not  at 
present  correctly  predict  many  of  the  radiative  and  dynamic  features 
observed  in  implosions.  Thus  the  first  theoretical  task  in  devising 
the  ideal  machine-load  configuration  is  to  determine  the  reasons  for 
the  existing  disagreement.  In  the  present  section  we  summarize  a 
new  set  of  simulations  which  extend  previous  modeling  by  including 
the  effects  of  (i)  a  short  circuit  in  the  pulse  power  feed,  and  (ii) 
a  phenomenological  anomalous  resistivity  in  the  pinched  plasma. 
This  study  is  part  of  an  ongoing  effort  to  uncover  those  physical 
processes  which  will  bring  the  simulation  results  more  into  line 
with  the  experimental  data.  The  present  results  for  the  pinch 
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radii,  total  and  K-shell  yields,  and  K-shell  pulse  widths  are 
encouraging  but  still  do  not  exactly  match  the  observations. 

The  present  simulations  were  done  for  aluminum  wire  array 

implosions  on  the  Double  Eagle  pulse  generator  of  Physics 
International.  This  choice  was  motivated  by  the  recent  set  of 
experiments  on  this  machine  and  the  accompanying  excellent  data  set 
kindly  provided  to  us  by  C.  Deeney  of  Physics  International.  For 
the  numerical  simulations  ve  used  the  ZPIMPO  radiation- 

magnetohydrodynamic  code.  This  code  has  been  fully  describe  in  an 
NRL  Memorandum  Report  #6448.  Briefly,  the  ZPIMPO  code  is  a  single 
zone  ("0-D")  cylindrically  symmetric  hydrodynamic  code,  i.e.,  the 
plasma  dynamics  and  state  is  modeled  as  a  shell  with  uniform  density 
and  temperature  and  separate  momenta  equations  for  the  inner  (R^) 
and  outer  (R  )  radii.  The  diffusion  of  the  magnetic  field, 
however,  is  followed  over  several  zones  in  order  to  model  the 
current  penetration  into  the  plasma  and  resistive  heating  correctly. 
The  ionization  state  of  the  plasma  in  the  shell  is  followed  by 
solving  time  dependent  rate  equations.  The  radiation  transport  uses 
the  probability  of  escape  formalism  for  the  lines  and  edges  and 
multi-frequency  for  the  free-free  emission.  The  sole  advantage  of 
this  single  zone  modeling  compared  to  a  full  multi-zone  1-D  code  is 
the  efficient  running  time  and  consequent  ability  to  survey  a  large 
range  of  parameter  space  in  a  cost  effective  manner. 

The  magnetohydrodynamics  in  the  plasma  is  coupled  self- 

consistently  to  an  equivalent  circuit  model  for  Double  Eagle  shown 

in  Fig.  1.  The  voltage  source  supplied  by  the  generator  VS0(jrce  is 

a  known  function  of  time,  while  the  feed  from  the  insulator  stacks 

to  the  load  is  characterized  by  fixed  inductances  La  and  L^,  and  a 

fixed  resistance  R  .  In  the  discussion  presented  below  we  have 

considered  the  effect  of  a  short  circuit  in  the  feed  with  a  variable 

resistance  R  .  The  location  of  the  short  divides  the  total  feed 
s 

inductance  into  a  part  before  the  short  and  a  part  after  the  short. 
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The  currents  to  be  determined  in  the  circuit  are  denoted  by  I  in 
Fig.  1.  The  imploding  plasma  provides  the  load  voltage  VjQacj. 

Ve  considered  three  sets  of  runs,  comprised  of  six  cases  each. 
The  first  set,  termed  the  standard  set,  has  no  short  in  the  feed  and 
uses  classical  Spitzer  resistivity.  For  the  second  set  ve  added  a 
short  in  the  feed  once  the  load  voltage  exceeded  3  MVolts.  The 
third  set  also  has  a  short  but  in  addition  ve  employed  a 
phenomenological  resistivity  vhich  increases  as  the  plasma  pinches. 
Specifically 


^Spi tzer x 


[“•[‘-(Wn]- 


R  .(0) 
out 


(*) 


The  parameters  for  the  three  different  run  sets  are  summarized  in 
table  I. 


TABLE  I 

SET  PARAMETERS 

run  set  short  resistivity 


A 

B 

C 


no  Spitzer 

yes  Spitzer 

yes  eqn.(*)  vith  3=100 
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Since  the  load  is  an  array  of  solid  vires  the  initial  conditions 
in  a  cylindrically  symmetric  model  are  the  total  mass  of  the  array 
and  its  radius  R(0).  The  initial  conditions  for  the  six  cases  in 
each  run  set  are  listed  in  table  II. 


TABLE  II 

INITIAL  CONDITIONS 


case 

Mass 

R(0) 

M*R(0) 

Ugm 

cm 

Ugm*cm 

1 

2222 

0.3 

200 

2 

556 

0.6 

200 

3 

200 

1.0 

200 

4 

89 

1.5 

200 

5 

50 

2.0 

200 

6 

22 

3.0 

200 

Figure  2a  presents  the  time  development  of  several  quantities 
for  case  3  of  the  run  set  A.  From  top  left  to  bottom  right  within 
each  figure,  we  show 

(i)  the  outer  (solid)  and  inner  (dashed)  radii; 

(ii)  the  feed  current  I  (dashed)  and  the  load  current  I,  (solid); 

(iii)  the  K-shell  (solid)  and  total  (dashed)  radiation  pulse; 

(iv)  the  electron  temperature  Tg; 

(v)  the  mean  charge  state  of  the  plasma  <Z>; 
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(vi)  the  total  (solid)  and  radiated  (dashed)  energy  in  the  pinch. 
Likewise,  Fig.  2b  presents  the  time  development  for  case  3  of  run 
set  B,  and  Fig.  2c  presents  the  time  development  for  case  3  of  run 
set  C. 

Figure  3a  presents  the  total  (filled  circles)  and  K-shell 
(circles)  yields  derived  from  the  numerical  simulations  for  the  six 
cases  of  table  II  in  run  set  A  as  a  function  of  the  initial  outer 
radius  R(0).  For  comparison,  ve  also  display  the  total  (filled 
triangles)  and  K-shell  (triangles)  yields  measured  on  the  Double 
Eagle  experiments.  Likewise,  Fig.  3b  presents  the  yields  from  run 
set  B,  and  Fig.  3c  presents  the  yields  from  run  set  C.  Again  the 
same  experimental  data  is  shown  in  each  figure. 

Ve  first  note  that  the  implosion  times  from  the  simulations 
range  from  85  to  95  nanoseconds  which  agree  with  the  experimental 
values  to  within  -  5  nanoseconds.  In  the  experiment  the  initial 
time  is  measured  by  extrapolating  the  linear  portion  of  the  current 
rise  back  to  I  =  0.  The  time  of  pinch  assembly  on  axis  is  denoted 
by  the  appearance  of  the  K-shell  pulse.  The  difference  between 
these  two  fiducial  times  is  the  implosion  time.  This  procedure  is 
indicated  in  Fig.  2a.  This  agreement  indicates  that  most  of  the 
mass  in  the  wire  arrays  must  participate  in  the  implosion,  otherwise 
the  experimental  implosion  time  would  be  significantly  shorter  than 
the  simulation  results. 

Figures  2a  and  3a  indicate  several  features  of  the  set  A 
simulation  runs  which  are  in  disagreement  with  the  experiment. 
First  from  Fig.  2a,  the  minimum  pinch  radius  of  the  simulation  is 
much  smaller  than  the  observed  pinch  radius,  which  is  about  0.1  cm. 
Second,  also  from  Fig.  2a,  the  K-shell  pulse  width  is  <  1 
nanosecond,  while  the  observed  pulse  width  is  on  the  order  of  10 
nanoseconds.  And  third  from  Fig.  3a,  the  total  and  K-shell  yields 
are  significantly  larger  than  the  observed  values.  Instead  of 
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collapsing  to  a  pinched  configuration  of  one  tenth  the  initial 
radius  for  ~  10  nanoseconds  as  in  the  experiments,  the  simulation 
models  execute  a  semi-periodic  bouncing  at  -  0.01  cm  with  K-shell 
radiation  emerging  at  least  during  the  first  bounce.  The  long  term 
state  is  a  -  hundred  nanosecond  long  Bennett  equilibrium  radiating 
in  the  L-shell  of  aluminum.  This  is  the  theoretical  radiative 
collapse  phase.  Ve  note  these  features  are  characteristic  of  1-D 
simulation  radiation-hydrodynamic  multi-zone  codes,  and  are  not  a 
peculiarity  of  our  computationaly  faster  one  zone  modeling. 

The  discrepancy  in  total  radiated  energy  suggests  that  more 
energy  is  coupled  to  the  pinch  load  in  the  simulations  than  in  the 
actual  experiment.  In  order  to  reduce  the  amount  of  coupled  energy 
ve  supposed  that  a  short  circuit  is  triggered  in  the  feed  once  the 
load  voltage  exceeds  a  certain  value.  Since  the  implosion  time 
occurs  before  the  time  of  peak  machine  current,  a  short  in  the  feed 
located  upstream  of  the  load  would  effectively  limit  the  energy 
reaching  the  load  by  reducing  the  load  current.  Run  set  B  contains 
the  simulations  with  a  short,  but  again  with  Spitzer  resistivity. 
The  results  are  displayed  in  Figs.  2b  and  3b.  From  Fig.  2b  ve  note 
that  the  current  upstream  of  the  short,  where  the  current  monitor  is 
placed  in  the  experiment,  continues  to  rise  after  implosion  of  the 
pinch.  The  current  in  the  load,  however,  drops  from  3  to  0.5  MAmps 
at  implosion.  The  total  and  K-shell  yields  of  the  set  B  runs  are, 
as  suspected,  smaller  than  the  yields  without  a  short.  (Compare  Fig. 
3b  for  set  B  with  Fig.  3a  for  set  A.)  As  a  matter  of  fact,  both  the 
total  and  K-shell  yields  are  on  the  order  of  the  experimental  values 
for  set  B  runs  with  a  short.  This  result  does  not  appear  to  be  very 
sensitive  to  changes  in  the  parameters  which  characterize  the  short, 
namely,  the  load  voltage  which  triggers  the  short,  the  splitting  of 
the  feed  inductance,  and  the  variation  of  the  short  resistance  Rg . 
The  existence  of  the  short  circuit  could  be  checked  by  using  several 
current  monitors  over  the  length  of  the  feed  and  looking  for  current 
profiles  like  those  of  Fig.  2b.  From  Fig.  2b  one  sees,  however,  that 
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the  pinch  radius  is  still  much  smaller  than  the  observed  0.1  -  0.2 
cm  and  the  K-shell  pulse  length  is  still  much  shorter  than  the 
observed  -  10  nanoseconds. 


Clearly  some  physical  mechanism  not  considered  in  the 
simulations  appears  to  cause  a  rapid  thermalization  and  a 
hydrodynamic  bounce  at  -  0.1  cm,  which  is  much  larger  than  the 
bounce  radius  found  in  the  simulations.  Vithin  the  context  of  a  1-D 
simulation  code  ve  can  mimic  the  effects  of  a  rapid  thermalization 
by  considering  a  phenomenologically  enhanced  resistivity  according 
to  eqn.(*).  Let  us  estimate  the  value  of  the  parameter  g  in  this 
model  equation  by  comparing  the  classical  resistive  heating  rate 
with  the  observed  mean  K-shell  power.  The  latter  quantity  is 


P  ,  (K-shell)  -  =  3xl01!  V 

rad  10  nsec 

From  the  observed  pinch  radius  (0.1  cm)  and  peak  current  (3  MA), 
assuming  the  Bennett  temperature  and  Spitzer  resistivity,  one  finds 
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Thus  the  classical  heating  rate  is  nearly  two  orders  of  magnitude 
smaller  than  the  K-shell  radiated  power.  Hence  we  considered  the 
effect  in  run  set  C  of  an  arbitrarily  enhanced  resistivity 
according  to  eqn.(*)  with  g  =  100.  The  results,  displayed  in  Figs. 
2c  and  3c,  do  succeed  in  limiting  the  pinch  radius  to  -  0.1  cm  at 
least  on  the  first  bounce,  producing  broader  K-shell  pulse  widths  ( 
-  5  nanoseconds)  compared  to  previous  runs,  and  predicting  total 
yields  in  better  agreement  with  the  experimental  values  than  the 
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previous  runs.  However,  the  K-shell  yields  have  increased  from  run 
set  B  to  about  three  times  the  measured  values  at  R(0)  =  1.0  cm. 

The  general  improvement  in  the  comparison  of  simulation  results 
with  experimental  values  provided  by  the  enhanced  resistivity  of 
eqn.(*)  is  encouraging  for  developing  a  predictive  capability. 
This  enhanced  resistivity  effectively  prevents  the  radiative 
collapse  observed  in  previous  simulations  by  continually  diffusing 
the  magnetic  field  into  the  plasma  and  increasing  the  resistive 
heating.  The  actual  physical  mechanism  behind  the  phenomenological 
resistivity  cannot  be  ascertained  from  the  present  1-D  model.  The 
enhancement  may  arise  from  convective  or  sausage  instabilities  in 
the  plasma  shell  during  the  run  in  phase.  The  form  of  the 
enhancement  factor  in  eqn.(*)  is  of  course  not  unique.  If  instead 
of  eqn.(*)  we  use  the  following  expression  for  the  enhanced 
resistivi ty 


"  *  Vtzer*  [  1  *  mi"(  e,<"sec)/25-  1  ,  100  )  ]  ,  (**) 

we  find  nearly  identical  results  to  those  of  run  set  C.  This 
suggests  that,  if  instabilities  are  responsible  for  the  enhanced 
resistivity,  their  growth  rate  should  be  on  the  order  of  25 
nanoseconds.  To  pursue  the  problem  of  understanding  the 
limitations  to  radiative  collapse  in  a  Z-pinch  one  must  investigate 
the  implosion  and  radiative  dynamics  with  a  multi-dimensional 
simulation  code. 
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Figure  3c 
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Exploiting  Radiative  Collapse  in  a  PRS 


Radiative  collapse  has  been  discussed  at  some  length 7-10  and,  while  the  phe¬ 
nomenon  has  not  been  observed  to  date,  the  concept  has  an  attractive  simplicity 
when  considered  alonside  the  complex  technological  requirements  of  imploding  PRS 
load  schemes.  A  PRS  made  to  exhibit  radiative  collapse  would  probably  consist  of 
several  parallel  cryogenic  fiber  sources,  made  either  of  pure,  moderate  Z  frozen 
gases  (Ne,Ar),  or  of  high  Z  seeded  hydrogen  fibers.  The  parallel  connection  would 
require  vacuum  convolutes  similar  to  those  contemplated  for  Reflex  Triode  bremm- 
strahlung  sources,  and  could  be  driven  directly  (on  ps  quarter  cycle  timescales) 
by  fast  capacitor  banks  like  those  used  in  the  Falcon  or  FastCap  decade  machine 
designs. 

In  view  of  the  conjectural  aspect  of  the  whole  "radiative  collapse ”  scenario, 
a  detailed  machine  or  even  front  end  design  would  be  premature,  but  one  can 
fruitfully  examine  the  source  dynamics  with  simple  models  in  order  to  crysfaiize 
more  fully  the  various  tradeoffs  involved.  IVe  will  do  this  here  for  two  cases.  The 
first  is  a  ” classical ”  radiative  collapse  at  the  Pease-Braginskii  current;  the  second, 
a  screwpinch  variant  on  the  scheme  which  allows  one  to  select  the  critical  current 
level,  or  at  least  vary  it  about  the  Pease-Braginskii  value. 

To  date  u-e  have  mapped  out  only  the  essential  features  of  these  schemes,  but, 
whenever  more  concrete  experimental  verification  of  the  phenomenon  is  available, 
the  model  will  readily  allow  investigations  into  the  mutual  interaction  of  such  paral¬ 
lel  loads  in  a  decade  machine  environment ,  or  the  expected  radiation  pulse  history, 
spectrum  and  yield. 
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Essential  features  of  a  classical  radiative  collapse 

The  refrigeration  is  driven  by  a  competition  between  ohmic  heating  and  radia¬ 
tive  cooling  in  a  Bennett  equilibrium.  For  such  an  equilibrium  pinch  large  current 
forces  a  density  large  enough  to  cool  the  plasma  quickly.  When  the  current  in  the 
pinch  exceeds  Ipb,  then 

A  _  ( Ph  ~  Pr)  _  .  Ipb_  >2  . 

a  ~  1 2  4/r  [  I  ~  \  ' 

predicts  a  rapid  contraction  of  the  pinch.  It  is  this  result  which  implies  the  use 
of  multiple ,  isolated  loads  in  parallel ,  since  the  pinch  current  will  not  significantly 
exceed  the  Pease-Braginskii  current  1.5  MA)  and  the  energy  sources  will  provide 
tens  of  megamperes. 

A  more  delicate  constraint  has  been  pointed  out  by  Robson11.  The  load  must 
accept  energy  for  a  sufficiently  long  time ,  viz.  a  collapse  time(tc )  =  several  radiation 
cooling  times  (tr).  A  density-  collapse  time  product  then  arises 

„  y/T 

n0tc  =  1.42  •  1013  - T -  . 

H  hereupon  the  use  of  a  high  initial  density  is  required ,  since  the  collapse  time  must 
be  short  enough  to  be  compatible  with  pulse  power  sources.  Here  the  temperature 
T  is  also  deremined  by  the  current  through  the  Bennett  condition. 

In  the  laboratory,  what  conditions  might  help  obtain  radiative  collapse?  Frozen 
Hn  fiber  pinches  are  observed  to  be  stable  so  long  as  the  Poynting  flux  is  inward. 
Typical  stable  domains  are  300  ns  to  0.5  MA.  An  "inductive  notch"  during  the 
compression  will  therefore  not  cause  disruption  until  the  voltage  has  reversed  after 
stagnation.  Since  the  pinch  and  its  immediate  inductance  are  a  final  energy  store 
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very  close  to  the  load,  the  use  of  relatively  slow  driver  timescales  (ft: l—*2  gs)  allows 
an  easier  energy  transfer  to  this  inductive  store.  The  spontaneous  collapse,  when 
it  occurs  (at  I  >  Ipb ),  would  allow  a  nearly  complete  conversion  of  this  nea rby 
magnetic  energy  to  radiation. 

What's  a  favorable  machine  configuration  for  obtaining  these  conditions?  A 
load  of  this  sort  would  need  no  opening  switch,  since  the  pinch  timescales  favor  a 
direct  drive.  The  quarter  cycle  time  (%1  ps)  for  a  Falcon  or  FastCap  bank  is  right 
in  a  useful  window  for  this  load. 

Parallel  loads  convoluted  at  the  usual  radius  of  the  ring  diode  in  bremsstrahlung 
mode  offer  the  best  energy  coupling.  The  overall  impedance  of  the  array  can  remain 
ft:  0.1  Oil  or  less  even  when  collapsed  to  very  small  radius.  Moreover,  the  individual 
load  modules  tend  to  be  self  correcting  in  collapse  time  because  increased  motional 
impedance  in  one  holds  the  current  down  while  others  catch  up.  These  considera¬ 
tions  lead  to  the  conceptual  design  shown  in  Figure  2. 

Fig.  2  A  Generic  Design 


Circuit  Element  Model  for  a  Bennett  PRS 

A  quasi-equilibrium  model  for  the  radial  motion  is  summarized  by 


a  ,  „ 
-  =  107ZU 


3/ 

21  ' 


Ipb 


=  433 


1  +  Z 
2  Z 
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In  our  generic  design  the  single  module  inductance  and  impedance  are  given 

by 

L  =  1.8  •  1021L/c2  In  —  [nH]  , 

a 

z  =  z0  -  -  1.8-10 12/p/c2  [n] , 

a 

and 

z0  =  i.io9.io%(A)^L  [n|. 

What  results  can  be  expected  in  such  a  radiative  collapse?  Consider  as  ” middle 
of  the  range ”  tes<  case  a  40  pm  diameter  fiber,  driven  by  a  1.1  pF  water  capacitor, 
with  a  600  ns  quarter  cycle  time.  The  initial  store  energy  is  309  kJ  and  the  initial 
pinch  line  density  is  7.54- 1017  ions/cm3 ,  the  resulting  evolution  is  shown  in  Figure 
3(a). 

Fig.  3a  Radiative  Collapse  of  a  Model  Fiber  Load 

The  radiated  energy  per  module  is  decomposed  into  three  classes;  we  find  the  yields 
to  be  (i)  total  radiation:  146  kJ,  (ii)  radiation  >5  keV:  26  kJ,  and  (Hi)  radiation 
>10  keV:  5  kJ,  respectively.  The  radiation  power  development  is  shown  in  Figure 
3(b).  Hence  between  ten  to  fifteen  such  modules  would  be  a  reasonable  array  from 
the  standpoint  of  power  flow. 

Fig.  3b  Radiated  Power  in  Three  Channels 


B.  Modifications  when  axial  fields  are  included 


T I ME :  1 00n  s/d i v 

Figure  3(a) 
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RADIATED  POWEF  JN  THREE  CHANNELS 


The  use  of  an  axial  field  component  can  perhaps  allow  some  control  over  the 
onset  of  this  spontaneous  collapse,  as  well  as  the  gross  stability  of  the  PRS.  An 
arbitrary  mix  of  exterior  and  interior  axial  field  strengths  is  allowed,  again  using 
the  self-  similar  model.  The  mix  can  be  adjusted  by  twisting  the  return  current 
path  to  Nwind  per  cm,  so  that 

Si  OC  (  f  N  u’ind  / I Bennett )  i  ^  2  COTlSt.  . 

Circuit  Element  Model  for  the  Screw  Pinch  PRS 
The  momentum  equation  is 


2*4  ,  l'k  vk 


U{o)  =  r0(— ^  +  -  3-  ~) 


a  a°  a 


with  a  =  r(r)/r0,  0,.  =  Te  4- ( f  )U,  so  that  Tt  =  0e/B(Qr)  and 


",  =  y/{Z  +  1)T(0e)/mtonrl  , 


cy/T 


I(t) 


al\/T^m~n0 


If  one  neglects  ohmic  heating,  radiation  and  ionization,  viz.  B(Qe  — >  1,  the  implied 
adiabatic  heating  would  replace  — 1 •  in  momentum  equation. 

The  energy  equation  becomes 


a 


Qe  =  -(4/3)T(ee)(-)  +  nji/oy  -  prad(o,0e) 

a 


with  a  radiation  model  given  by 


P re.  ,v>,0f)  =  Min 


Pcre(Oc  ) 


a‘ 


qTA(Qr) 

1  +  2(0,) 
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here  the  "grayness”  is  modeled  by  G(Oe)  =  Pc  PE  /  Pbb  &t  the  constant  density  used 
to  calculate  the  CRE  emission  power  and  equation  of  state.  The  inductance  and 
impedance  now  include  the  return  current  windings  explicitly 

L  =  1.8  ■  10 21/p/c2  (In  —  +  I(  r2  _  r2 ) [ 2 7T Ar ] 2 )  [nH\  , 

r  Z 

Z  =  np  -  (  -  +  \2nN}2fr)  1.8  •  1012ip/c2  [«]  , 
r 

np  =  3.257  .l0'6(lp/nrl)(^  )-{-p^~  [Q]  . 

10  a277(©e) 

When  the  implosion  is  adiabatic  an  effective  potential  determines  the  turning 
points 

V(a)  =  4,(///b)  +  + 

The  kinetic  energy  of  the  pinch  is  then 

T2(a)  =  T;  -  V(a)  . 

Because  of  the  1/a3  scaling  of  axial  field  stress  the  interior  turning  points  persist 
until  larger  I / Ig,  Aru,mrf>  or  T2  values  can  erase  the  dip  in  T(o),  as  shown  in  Figures 
4(a,b,c ).  The  extinction  of  the  adiabatic  turning  points  with  larger  A  u>md  is  shown 
in  Figure  5.  A  critical  value  of  Nu.,nd  0.08  is  obtained  beyond  which  no  rebound 
is  possible  -  the  load  is  forced  to  a  singular  compression. 

When  Ohmic  heating  is  included ,  along  with  ’’gray”  radiation,  the  critical  value 
of  Nunnd  increases  substantially,  but  the  hole  is  still  there.  A  scan  of  Nwinj  from 
0.0— >1  is  enough  to  show  the  effect.  For  Nu.ind  ~  1,  a  is  ever  more  negative  for 
a  <  0.749,  even  though  ©<•  remains  positive  as  the  load  collapses.  For  Nwt„d  ~  0,  a 
is  ever  more  positive  and  the  first  turning  point  occurs  for  a  =  0.727.  On  or  above 
Nunnd  ~  0.8,  a  radiative  collapse  is  unavoidable. 
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Figure  4a 


Adiabatic  Potential  for  Initial  Kinetic  Energy.  to= 1 .0  .  Nu.ind  =  0.125 
Uariabie  l/lb  =  [1.0.  6.0] 
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Figure  4b 


Adiabatic  Potential  for  Initial  Kinetic  Energy.  to=.10  .  I/lb  =  1.0 
Uariable  Nujnd=]0.  0 .  0.25] 
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Figure  4c 


Adiabatic  Potential  for  N urn d =0.1 25.  I/lb  =  1.0 
Uariable  Initial  Kinetic  Energy.  to=[0 .0 . 2.5] 
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ING  POINT 


TURNING  POINTS  AS  A  FUNCTION  OF  WINDING  NUMBER 

PER  UNIT  LENGTH 


Conclusions 


Many  inductive  store  concepts  separate  fuse ,  switch  and  load;  others  make  the 
load  and  switch  one  unit.  In  the  sense  discussed  here,  because  of  the  critical  trigger 
current  for  collapse,  the  fuse  can  be  made  the  load.  The  Pease-Braginskii  current 
is  only  one  limit  on  the  current  per  load  module  and  this  natural  trigger  value  for 
highly  efficient  energy  conversion  is  perhaps  tunable  with  axial  field  components. 
Such  a  PRS  device  allows  one  to  charge  the  final  inductive  store  on  bank  quarter 
cycle  timescales  and  recover  the  energy  quickly  in  only  one  stage  of  power  ampli¬ 
fication.  A  direct  drive  and  parallel  load  modules  offer  a  ready  compatibility  and 
a  straightforward  development  path  with  presently  conceived  decade  machine  de¬ 
signs.  High  energy  (hu  >  10  keV)  yields  on  the  order  of  5  kJ  per  load  module 
could  be  realized. 
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IV.  PRS  Optimization:  Electron  and  Field  Dynamics  in  the  z-pinch 

Introduction  and  Review 

Energetic  electrons  can  be  very  influential  in  the  dynamics  of  a  z-pinch  implosion  as 
well  as  in  the  radiation  generated  by  the  pinch,  yet  this  influence  has  not  been  investigated 
in  any  depth  to  date.  Distortions  of  the  electron  distribution  function  from  a  Maxwellian 
shape  will  change  transport  coefficients  (e.g.,  conductivity)  and  alter  rates  of  radiation 
mediated  by  electron-ion  collisions.  Beams  of  electron  runaways  such  as  have  been  observed 
seem  to  be  connected  with  the  “bright  spots”  from  which  a  large  part  of  the  pinch  K-shell 
radiation  typically  emerges.  Bremsstrahlung  from  energetic  electron  beams  can  also  be  a 
significant  source  of  hard  x-rays.  With  increased  machine  power,  these  effects  could  be 
even  more  significant  and  their  understanding  more  important. 

The  goal  of  this  study  is  to  understand  the  role  of  energetic  electrons  in  a  z-pinch 
implosion.  In  earlier  work,  the  Fokker-Planck  equation  was  analyzed  to  determine  factors 
significantly  affecting  the  electron  distribution  function.  A  one- dimensional  numerical 
solution  of  the  Fokker-Planck  equation  was  developed  and  used  as  a  diagnostic  of  one¬ 
dimensional  MHD  z-pinch  simulations.  From  this,  estimates  of  the  high-energy  electron 
population  and  degree  of  anisotropy  were  obtained  at  all  times  in  the  simulated  implosions. 

Three  factors  affecting  the  electron  distribution  function  were  studied:  ohmic  heating, 
charge-conserving  inelatic  collisions,  and  compression.  Ohmic  heating  was  found  to  be  the 
most  important  in  changing  the  shape  of  the  distribution.  The  inelastic  collisions  studied 
were  generally  negligible,  and  compressional  heating  was  largely  isotropic,  even  when  of 
significant  magnitude  (as  at  stagnation),  and  had  the  principal  effect  of  accentuating 
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ohmically-produced  anisotropies. 

Ohmic  heating  causes  an  enhancement  of  the  distribution  function  focused  at  the 
electron  energy 


where  ve{  is  the  thermal  electron-ion  collision  frequency,  and  ujc  is  the  electron  cyclotron 
frequency.  The  magnetic  field  determines  the  character  of  the  heating:  for  low  magnetic 
fields  lu'c  <  2uei,  the  superthermal  tail  is  enhanced,  while  for  higher  fields,  the  tail  is 
depleted.  The  magnitude  of  heating  is  determined  by  the  electric  field:  with  a  zero 
magnetic  field,  significant  numbers  of  runaways  are  generated  by  the  Dreicer  field 

Ep  =  -vte\/Z  mkT, 
e 

where  e  and  m  are  the  electron  charge  and  mass,  and  uee  is  the  thermal  electron-electron 
collision  frequency.  With  a  magnetic  field,  a  higher  electric  field  is  needed  to  produce 
energetic  electrons,  and  the  electrons  generated  are  more  isotropic  and  less  like  classical 
“runaways”.  A  strong  magnetic  field  (as  defined  above)  effectively  shuts  off  runaway 
production  and  focuses  the  system’s  energy  on  the  well-confined  bulk  of  the  distribution. 

A  system  that  generates  a  significant  number  of  energetic  electrons,  then,  must 
experience  a  large  electric  field  E  >  Ep  while  the  magnetic  field  is  still  relatively  small 
(u/c  <  2uet)-  A  burst  of  directional  electrons  should  result,  which  can  be  amplified  by 
compressional  heating,  should  it  occur  at  that  time.  Energetic  electrons  are  rapidly 
confined  by  a  large  magnetic  field,  and  high  electric  fields  cannot  overcome  this  effect. 

Only  a  small  runaway  electron  population  was  predicted  from  the  z-pinch  implosion 
model  used  for  diagnostics.  In  particular,  the  electron  population  at  stagnation,  while  non- 
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Maxwellian,  was  found  to  be  highly  isotropic.  This  is  not  in  accord  with  experiment,  where 
bursts  of  electrons  have  been  seen  at  stagnation.  The  reason  for  the  discrepancy  is  believed 
to  be  incomplete  knowledge  of  the  radial  distribution  of  the  electric  and  magnetic  field 
strengths  inside  the  pinch.  Small  changes  in  the  model  field  distributions  could  significantly 
affect  the  predicted  runaway  populations.  There  is  little  experimental  information  about 
pinch  field  strengths,  which  depend  strongly  on  the  way  the  pinch  current  is  carried.  A 
number  of  different  theoretical  models  of  pinch  current  exist,  and  there  have  been  some 
interesting  numerical  results  in  the  related  topic  of  plasma  opening  switches  that  suggest 
that  magnetic  field  diffusion  is  much  faster  than  a  Spitzer-based  fluid  theory  would  predict. 
With  plasma  opening  switches,  initial  conditions  are  very  important  in  determining  the 
evolution  of  the  system. 

Plasma  opening  switches  utilize  much  lower  densities  than  the  z-pinch,  especially  from 
an  exploding  wire.  The  parametrizing  study  of  the  pinch  corona  presented  here  is  a  first 
step  in  studying  the  factors  determining  current  and  field  distributions  in  an  imploding 
z-pinch. 

The  coronal  current  channel  in  an  exploding-wire  z-pinch 

The  corona  formed  outside  an  exploding  wire  plays,  at  least  initially,  an  important 
role  in  current  conduction  and  pinch  dynamics.  Because  of  its  high  conductivity,  the 
pinch  medium  tends  to  exclude  magnetic  field  and  current.  In  the  corona,  however, 
turbulence  will  limit  the  conductivity  and  an  anomolous  resistivity  will  allow  rapid  field  and 
current  penetration.  The  low  density  of  the  corona  limits  radiative  losses  and  allows  high 
temperatures  to  be  attained  which  paradoxically  enhance  the  conductivity.  So  the  corona 
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can  carry  a  large  fraction  of  the  pinch  current  while  also  maintaining  a  quite  different 
configuration  of  electric  and  magnetic  fields  from  the  pinch  core. 

In  addition  to  its  implications  for  electron  energy  transport  and  radiation,  the  corona 
potentially  has  a  significant  effect  on  the  overally  pinch  dynamics,  since  it  acts  as  a 
competing  current  channel  to  the  core  and  heats  the  core  by  energetic  electron  collisions. 
These  effects  are  incorporated  into  most  pinch  models  by  adding  a  thin  skin  current  layer 
to  the  current  channel. 

To  study  the  behavior  of  the  corona  under  various  pinch  conditions,  parameters  from 
a  ID  MHD  z-pinch  model  were  used  to  determine  the  dynamical  evolution  of  the  corona 
extent,  density  and  temperature.  The  parameters  used  were  appropriate  for  a  titanium  wire 
implosion.  The  simplest  case  is  examined:  a  corona  where  limited  conductivity  permits 
field  penetration  and  large  currents.  Low  densities  minimize  radiative  losses,  and  very  high 
temperatures  and  currents  result.  Initial  theoretical  work  on  this  problem  was  performed 
by  D.A.  Tidman. 

Z-pinch  model: 

The  system  is  considered  to  be  cylindrically  symmetric  with  spatial  uniformity  in  the 
z  direction.  The  ion  atomic  mass  number  is  A ,  and  degree  of  ionization  is  Z.  Corona 
parameters  are: 

a  =  Corona  inner  radius/core  outer  radius 

7j  =  Core  current  (amp)  72  =  Corona  current  (amp) 

kT\  =  Core  temperature  (eV)  fcT2  =  Corona  temperature  (eV) 
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The  corona  is  distinguished  from  the  core  by  having  a  limited  conductivity,  such  that: 


j  =  ec3ne 

where  nt  is  the  electron  number  density,  and  ca  is  the  ion  acoustic  sound  speed: 

ca  =  9.788  x  10s  \JZkTi) A  <  vth 

The  core  dynamics  are  determined  by  MHD  model  (with  Spitzer  resistivity  and  line 
radiation).  Pinch  is  divided  into  30  uniform  annular  cells,  numbered  from  the  inside 
out.  Outside  the  outermost  cell  is  the  isothermal  corona. 


Fluid  Equations: 


The  complete  fluid  equations  are  listed  here  for  reference.  These  are  approximated 
for  the  corona  considered  here  by  neglecting  the  right-hand  side  of  each  equation.  The 
assumption  is  that  the  corona  is  quasi-static,  following  the  motion  of  the  core,  and  that 
it  is  in  Bennet  equilibrium.  This  is  just  the  simplest  assumption  for  a  first  look  at  the 
problem. 


1-  _  dV  /-  ~\  -  -  dV 

-j  xB-V  P  =  p~  +  p(V-V)V  +  neE- p— 


coll 


-  -  47r  -  1  dE 

VxB - j  = 

c  c  crt 


dt  \2 

Quasi-Static  Corona: 


d  (3  ,  \  -  *  ^  ^  d  (\ 

-(-nkTj=Ej-Vq-m^- 


nmV 2  ]  +  —  ( —  nrrtl'2 

dt  V  2 


coll 


The  equations  actually  satisfied  by  the  corona  are: 

--jB  =  kT2 ^ 
c  2  dr 

c  1  d 

T*-rTr(rB)=3 
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These  can  be  solved  analytically.  The  marginal-current  approximation  relates  j  and  n,  so 
that  one  equation  for  (ri?)  can  be  obtained: 


J_  ]_d_  ( r2B 2 
47r  r2  dr 


ckT2  d 
4neca  dr 


7^(rB) 


The  first  integration  of  this  expression  (from  a  to  r)  gives  an  equation  for  y  =  rB: 


dy  ~ 

r—  =2 y  -  ay2  +  /3 


The  second  integration  gives  an  expression  for  y(r): 


y  = 


V2  +  (r/ro)26yi 
1  + (r/r0)26 


where 


Vi,2  =  —(1  ±  \/l  -  fia) 

a 

6  =  a/1  +  0a 


This  expression  for  y(r)  is  plotted  in  fig.  (1),  and  the  corresponding  magnetic  field  is  shown 
in  fig.  (2).  These  figures  show  the  shape  of  the  quantities  only,  which  remains  the  same 
for  any  values  of  the  parameters.  Different  parameter  values  will,  of  course,  change  the 
scaling. 
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Profiles: 


The  scale  Length  r0  is  given  by 


The  density  profile  n(r)  is  then: 


T  +  ( a/r0 )26'  2  (r/r0)26  2 
1  +  (r/r0)26_  ( a/r0)2i~ 2 


The  density  profile  n(r)  is  plotted  in  fig.  (3).  The  scale  will  change  but  the  shape  will 
remain  invariant  for  different  parameter  values. 


Heating 


term 


To  get  an  accurate  estimate  of  the  corona  conductivity  and  energy  interaction  with 
the  core,  it  is  important  to  have  an  accurate  heat  transport  equation.  This  is  one  of  the 
most  difficult  problems,  since  large  amounts  of  heat  are  exchanged  across  the  core-corona 
boundary,  and  since  heat  transport  across  large  thermal  gradients  such  as  axe  set  up  across 
the  boundary  is  as  yet  poorly  understood.  For  the  present,  the  simplest  approximation  is 
used,  namely  that  of  no  net  heat  transport  across  the  boundary.  The  is  certainly  true  at 
early  time,  but  thermal  conduction  from  the  hot  corona  to  the  core  is  likely  to  be  important 
and  influential  to  both  core  and  corona  at  later  times,  and  should  be  stud^d  in  the  future. 
The  coronal  ohmic  heating  term  is 

I  (jntr2)  =  J(r)  ■  E 


which  becomes,  since  j  =  encs : 


dkT2 

dt 
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Boundary  conditions: 

At  t  =  0,  T\  =  T2  and  n2(a)  =  nj(a).  The  electric  field  E  is  considered  to  be 
continuous  across  the  core-corona  boundary.  For  t  >  0,  current  is  continuous  across  the 
corona  boundary': 

Ji 

n0 - T^TT 

ecs(r2) 

yo  =  —  =  0.41!  [A] 
c 

y'o  =  a 

Note  that  yi  >  yo  >  2/2- 
Results 

Corona  characteristics  for  a  Titanium  z-pinch  implosion  are  shown  in  fig.  ( 4)— ( 11)  as 
a  function  of  implosion  time.  Parameters  from  the  MHD  simulation  are  first  shown:  The 
average  ionization  Z  and  the  corona  inner  radius  a  are  shown  in  fig.  (4);  the  electric  and 
magnetic  fields  are  shown  in  fig.  (5);  and  the  temperature  and  current  are  shown  in  fig.  (6). 
Then  the  corona  parameters  axe  shown:  Fig.  (7)  shows  the  scale  length  r0  compared  to  the 
corona  inner  radius;  fig.  ^8)  shows  the  theoretical  position  of  maximum  electron  density  (at 
a  radius  smaller  than  the  inner  corona  edge  a),  while  fig.  (9)  shows  the  density  at  the  inner 
edge  and  at  the  maximum;  fig.  (10)  shows  the  core  and  coronal  current,  and  fig.  (11)  shows 
the  core  and  coronal  temperature.  Note  that  the  coronal  current  is  nearly  always  at  least 
an  order  of  magnitude  greater  than  the  core  current.  The  later-time  conductivity  of  the 
corona  is  not  realistic,  however,  because  of  the  impossibly  high  temperatures  it  attains. 
This  is  because  there  is  no  cooling  mechanism  for  the  corona,  and  the  conduction  that 
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certainly  occurs  in  reality  is  not  accounted  for.  Thus,  only  the  early-time  results  should 
be  considered  for  dynamical  conclusions.  Still,  these  results  do  show  the  importance  of 
heat  conduction  in  cooling  the  corona,  and  motivate  further  work  on  finding  a  practical 
heat-transport  term. 

Conclusions 

•  The  structure  and  characteristic  parameters  of  a  current-limited  region  have  been 
presented.  In  this  model  corona,  current  is  limited  by  scattering  from  ion-acoustic 
instabilities. 

•  Time  evolution  was  followed  of  a  turbulent  corona  of  a  simulated  titanium  wire 
z-pinch.  Scale  length,  temperature  and  total  current  were  calculated  and  compared  to  the 
pinch  values. 

•  Although  conductivity  is  limited,  increased  penetration  of  the  electric  field  and 
reduced  magnetic  field  results  in  significant  heating,  which  enhances  conduction. 

•  The  corona  can  carry  a  significant  fraction  of  the  pinch  current,  especially  early  in 
the  pinch. 

•  Accurate  heat  transport  is  very  important  to  the  model.  Strong  heating  can  cause 
the  model  to  break  down.  Future  work  should  include  better  calculation  of  radial  heat 
conduction  into  the  core  as  well  as  losses  from  atomic  line  radiation. 
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V.  New  PRS  Approaches 


a)  Modeling  The  Fast  Screw  Pinch  PRS  Load 


Introduction 


Plasma  Radiation  Sources  (PRS)  have  always  presented  something  of  a 
dilemma  to  theorist  and  experimenter  alike  in  that  instability  channels  seem  to 
occur  more  frequently  than  not  in  experiments,  making  it  difficult  to  obtain  de¬ 
tailed  and  meaningful  comparisons  between  theory  and  experiment. 

Yet,  in  order  to  scale  K-sliell  yields  to  higher  values  in  moderate  Z  elements,  it 
is  highly  desireahle  that  stable  array  implosions  be  achieved  with  large  initial  array 
diameters,  which  then  allow  a  longer  rundown  period  for  instability  growth.  The 
known  radiative  efficiency  of  one  dimensional  PRS  implosions  in  converting  kinetic 
energy  to  heat  and  ( K-shell )  X-radiation  leads  therefore  to  efforts  at  stabilizing  the 
laboratory  plasma  and  thereby  achieving  the  theoretical  K-shell  yields  characteristic 
of  large  diameter  implosions.  There  is  no  shortage  of  data  indicating  that  both 
m  =  0  and  m  =  1  motions  are  present  in  PRS  loads,  but  the  most  disruptive  of 
these  is  the  "sausage"  mode.  The  general  strategy  in  sausage  stabilization  is  to 
provide  an  axial  magnetic  held  component  which  cannot  be  ” sidestepped ”  by  the 
plasma  through  an  interchange  of  field  and  fluid,  nor  can  the  axial  field  be  as  easily 
compressed  -  it  builds  ” backpressure ”  with  an  r~3  scaling  as  compared  with  the 
r-1  scaling  of  the  pure  z-pinch  stress.  The  next  order  of  refinement  of  this  basic 
idea  involves  the  details  of  the  magnetic  shear  profile  in  the  load,  for  in  the  simplest 
limit  it  is  not  so  much  the  absolute  magnitude  of  the  axial  field  as  the  rapidity 
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of  spatial  change  in  the  ’’safety  factor ”  wJiich  stabilizes  things.  If  one  can  find  a 
practical,  effective  balance  of  axial  and  azimuthal  magnetic  fields,  there  will  result  a 
close  laboratory  approximation  to  the  idealized,  one-dimensional  implosion  process 
and  a  wider  variety  of  load  options,  probably  more  forgiving  in  design  tolerances, 
becomes  accessible.  For  the  screw  pinch  PRS  then,  our  primary  challenge  is  finding 
a  means  of  suppressing  the  load  instabilities  in  spite  of  all  natural  tendencies  to 
the  contrary.  The  physical  model  is  more  conventional  and  homogeneous  in  its 
composition,  but  the  risk  is  that  such  an  implosion  path  may  not  exist. 


On  the  other  hand,  should  the  proclivity  of  Z  pinches  to  distort  into  ” sausages ” 
and  "kinks”  be  suppressed?  Several  approaches  can,  in  principle,  lead  to  the  ex¬ 
ploitation  of  these  instabilities  to  enhance  the  radiative  yield  -  perhaps  through 
non-thermal  processes1,2 .  A  successful  development  along  these  lines,  e.g.  using 
X-pincltes,  must  deal  with  very  delicately  balanced  systems,  designed  to  close  tol¬ 
erances  and  likely  to  develop  large,  rapidly  evolving  gradients  in  plasma  conditions. 
Because  of  these  stiff  gradients,  many  competing  processes  must  be  folded  into  the 
analysis.  X-pinch  loads  which  might  be  productive  sources  of  runaway  electrons 
evolve  into  twin  axial  blobs  with  the  current  sheet  expanded  in  a  more  diffuse 
structure  near  the  midplane.  A  hot  cavitated  region  between  the  blobs  would  be 
the  expected  source  for  useful  fluences  of  supra-tliermal  electrons.  Such  a  region 
would  be  only  weakly  magnetized,  semi-collisional  and  quite  likely  microturbulent 
as  well.  The  adjoining  denser  portions  will  be  magnetized  and  more  fluid  like  and 
thus  domains  requiring  quite  different  theoretical  limits  can  be  expected  to  be  in 
intimate  contact.  The  initial  conditions  required  to  achieve  the  bifurcated  end  state 
and  produce  a  reasonably  deep  density  cavity  are  based  on  specific  mass  loadings 
as  functions  of  radial  and  axial  position.  While  the  minimum  tolerances  in  this 
mass  arrangement  required  for  a  given  level  of  performance  are  not  known,  it  is 


101 


reasonable  to  expect  a  system  somewhat  unforgiving  of  perturbations  away  from 
the  preferred  initial  state.  It  is  a  fair  summary  that  all  phases  of  X-pinch  evolution 
contain  rather  challenging  physics  to  unravel  before  a  proper  design  is  constructed. 
If  all  the  details  can  be  properly  investigated  and  traded  off,  there  will  result  a 
kind  of  ” plasma  judo ”  in  which  the  natural  motions  of  the  load  are  amplified  and 
channeled  to  the  preferred  evolution ,  but  exploiting  and  fine  tuning  the  instability 
is  risky  because  so  much  disparate  physics  must  be  included  in  the  design  model. 

Here  we  will  explore  the  first  of  these  equally  promising  general  design  path¬ 
ways,  a  mix  of  magnetic  field  components.  The  blend  of  z  and  9  pinch  field  sets 
(respectively,  the  TM  and  TE  modes)  into  a  helical  magnetic  field  lias  been  explored 
extensively  in  the  CTR  community.  For  fusion  devices  global  stability  is  essentially 
the  only  issue  and  perturbations  must  be  damped  on  very  long  timescales  com¬ 
pared  to  th  PRS  problem.  In  contrast,  for  a  modern  pulseline  driver,  the  stability 
is  needed  for  100  ns  timescales  only  and  the  load  can  be  annular  for  much  of  the 
rundown.  Since  the  annular  load  has  the  advantage  of  avoiding  any  substantial 
backpressure,  it  is  (at  least  marginally)  stable  until  the  stagnation  and  compression 
phase  wherein  the  timescales  are  futher  reduced  to  10  ns.  Since  the  energy  coupling 
efficiency  between  pulseline  and  load  is  improved  with  a  greater  dynamic  range  in 
load  inductance  during  the  rundown,  the  ability  to  achieve  a  stable  compression  over 
several  centimeters  of  radius  using  a  screw  pinch  load  would  be  extremely  beneficial 
to  the  PRS  effort  in  general. 

The  screw  pinch  could  therefore  be  an  attractive  PRS  load  if  (i)  it  achieves 
a  reasonable  measure  of  stability  against  m=0  deformations,  (ii)  it  can  gain  any 
advantage  in  power  coupling  to  modern  generators,  and  (Hi)  these  things  can  be 
achieved  with  little  or  no  loss  in  radiative  efficiency.  Using  a  radiatively  coupled, 
1-d  ideal  MHD  fluid  model  equipped  with  an  appropriate  equation  of  state,  the 


102 


global  stability  of  various  initial  shear  profiles  can  then  be  explored  using  modern 
criteria. 

Until  recently ,  the  local  stability  criterion  due  to  Suydam,  as  generalized  to  a 
screw  pinch  by  Hameiri3 ,  was  the  most  general  applicable  test.  It  is  now  established 
that  global  modes  can  exist  in  a  screw  pinch  under  some  conditions ,  and  the  as¬ 
sessment  of  stability  on  any  1-D  set  of  plasma  profiles  now  rests  on  several  distinct 
tests,  following  Bondeson4,5 , et.al. 

There  is  one  local  criterion  (generalizing  Hameiri's  result),  in  either  the  mag- 
netohydrodynamic  (MHD)  or  guiding  center  plasma  (GCP)  limit.  Bondeson  finds 
in  the  absence  of  axial  or  azimuthal  flow  that  the  ” modified  Suydam ”  spectrum  is 
stable  unless  at  arizy  r0 


Such  a  local  criterion  is  not  sufficient  for  stablility  in  the  ideal  MHD  limit  for  there 
global  modes  can  arise.  These  modes  are  stabilized  in  the  GCP  limit  but.  if  one 
has  an  MHD  profile  such  that  S(r0 )  =  (B2  —  yP)pih2  —  7P(k  •  B)2  =  0  an/yirikere, 
then  the  ”, global  slow  modes ”  are  unstable  if 

Ds(—,S",B,B,1k-B,k)>]  . 
vA  4 


Alternately,  if  one  has  a  MHD  profile  such  that  A(r0)  =  pui2  —  (k  •  B)2  =  0  tvrvy- 
•uukere,  then  the  ”, global  Alfven  modes’’  are  unstable  if 

D A{  — ,  A" ,  B ,  B'z,  B'e,k  •  B)  >  \  . 

vA  4 

In  these  expressions  u>  is  the  local  Doppler  shifted  frequency,  u>  =  u>  +  mv$/r  —  kv~ ,  7 
is  the  adiabatic  index,  while  p  is  the  winding  number  (it ’s  derivative  is  the  magnetic 
shear).  It  is  related  to  the  usual  ” safety  factor”  by  qR  =  p~l  =  Bg/rB..  The 
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expressions  Da  and  Ds  are  complicated  functions  of  the  arguments  listed,  which 
we  do  not  need  to  examine  here,  and  for  simplicity  the  modifications  in  the  Suydam 
index  due  to  axial  and  azimuthal  flow  are  suppressed.  Those  modifications  are  quite 
important  in  the  MHD  limit,  however,  and  can  be  included  readily  as  the  model 
develops. 

From  this  discussion  it  is  worth  noting  that,  in  the  MHD  limit,  we  must  deal 
with  spatial  derivatives  up  through  second  order  in  assessing  the  stability.  For 
this  reason  some  care  has  to  be  taken  in  building  an  accurate  solution  of  the  fluid 
equations  so  that  noisy  fluid  profiles  do  not  lead  to  erroneous  judgements  on  the 
worth  or  accessibility  of  any  particular  load  design.  A  present  only  one  screw  pinch 
load  experiment  is  active,  at  Maxwell  Laboratories,  Inc.;  we  will  examine  it  in  some 
detail  later.  The  range  of  possible  experiments  is  quite  wide,  spanning  all  manner 
of  axial  field  profiles,  and  is  largely  unexplored.  The  mapping  our  present  model 
makes  possible  is  intended  to  help  narrow  the  field  of  approaches,  but  to  do  that 
one  must  be  able  to  quantify  the  global  accuracy  of  the  fluid  solution  in  some  detail 
for  various  problems  and  timescales.  It  is  for  this  reason  that  the  benchmarking 
process  has  been  pursued  with  tenacity. 

We  will  present  below  some  variations  on  well-known  self-similar  MHD  flows 
which  admit  a  very  straightforward  (0-d)  "gas  bag"  model  that  in  turn  allows  us 
to  simulate  existing  screw  pinch  load  schemes  in  the  context  of  the  DNA  decade 
machine.  We  will  discuss  the  details  of  the  "gas  bag"  model  and  examine  the  energy 
coupling  tradeoffs  with  increasing  external  axial  magnetic  field  in  the  load  region. 
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The  Simplest  Picture  of  Energy  Coupling 


A.  The  Gas  Bag  Model 

The  self-similar  implosions  can  be  used  as  benchmarks,  but  they  also  have  a 
second  application  as  generic  pinch  models.  The  utility  of  such  models  is  clearly  lim¬ 
ited  to  the  identification  of  trade-offs  in  load  mass,  initial  position  and  the  like;  the 
caveats  associated  with  their  application  are  numerous.  In  PRS  theory,  whenever 
such  a  model  is  used,  we  are  implicitly  saying  that  either  (i)  the  detailed  profiles  of 
the  0-d  model  are  really  accessible,  or  (ii)  some  interior  dissipative  process  acts  to 
make  such  profiles  emerge,  or  (Hi)  they  represent,  some  radially  averaged  load  state. 
Here  a  mix  of  all  these  limits  is  contained  in  the  model. 

The  profiles  which  arise  are  expressed  most  conveniently  in  terms  of  the 
dimensionles  ratio  £  =  r/a  =  r0/a0;  homogeneous  compression  requires  that 
r  —  a(t)r0,and  a  =  a(t)a0.  Once  given  the  time  development  ofa(t)  the  profiles  are 
then  seen  to  ” breath ”  in  and  out  in  the  laboratory  frame  while  they  remain  frozen 
to  the  £  value  in  the  locally  advecting  frame  defined  by  t*r(£).  For  the  screwpinch, 
neglecting  normalizing  factors,  one  has 

n(£)  =  exp  - 

£>-(£)  =  ~  €XP  —  £2)  +  £2 exp  —  £2 

i(0  =  We  =  v/(l  -  (1  +  i2)cxp  -  i2) 

VA()  =  c*tr0 

Are  these  special  profiles  in  any  sense  accessible  to  the  load?  Plasma  viscos¬ 
ity  could  play  a  role  in  producing  and  maintaining  self-similar  states  if  it  is  both 
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nonlocal  (characterized  by  a  mean  relaxation  time  t)  and  much  stronger  than  clas¬ 
sical  small  angle  collisions  would  predict.  Clearly  the  current  programming  path 
to  homogeneous  compressions  is  not  genera]  enough  to  be  a  dominant  process,  so 
some  dissipative  mechanism  would  be  required.  A  second  ingredient  of  homoge¬ 
neous  compressions  is  an  isothermal  temperature  profile.  In  the  loosely  defined 
’’gray  body ”  limit,  viz.  sufficient  opacity  to  make  the  radiative  cooling  profile  peak 
near  the  surface,  the  interior  radiative  exchange  of  energy  is  usually  very  effective 
at  keeping  the  temperature  constant  over  the  bulk  of  the  pinch.  This  implies  that 
ohmic  heating,  if  included,  can  be  done  in  a  spatially  averaged  sense  without  accon- 
ting  for  the  detailed  (and  nonuniform)  current  profiles  implied  by  the  self-similar 
motion  -  only  form  factors  of  0(1 )  would  modify  the  resistance  in  that  case.  The 
Gaussian  density  profile  is  surely  not  the  only  one  likely  to  emerge  from  the  closure 
of  an  annular  PRS  load,  but  it  does  resemble  the  density  profiles  obtained  in  1-d 
calculations. 

On  the  other  hand,  the  resulting  ODE  for  U  =  a  possesses  the  proper  scal¬ 
ings  with  radius  (a)  for  all  the  separate  components  of  the  stress  -  pressure! v2s ), 
axial  magnetic  field  pushing  in  or  out  (v%t  ),  azimuthal  magnetic  field  ( i/gw ),  and 
viscosity(r).  If  instead  of  separating  the  momentum  equation  and  solving  for  pro¬ 
files,  we  had  spatially  integrated  the  separate  stress  components  and  used  the  av¬ 
erage  values  to  assign  a  strength  to  each  of  the  separation  constants,  the  result 
would  have  been  the  same.  In  other  words,  profiles  in  r0  which  differ  from  the  spe¬ 
cial  self-similar  ones  would  simply  produce  slightly  different  values  of  the  (arbitrary) 
separation  constants  if  averaged  over  some  spatial  domain.  The  details  of  such  form 
factors  are  hopefully  less  important  in  the  pinch  dynamics  than  the  radial  scalings 
of  each  competing  stress  component,  although  certainly  these  details  are  not  trivial. 
Only  comparisons  with  full  1-d  calculations  can  resolve  the  most  salient  of  these 
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differences  but  in  the  spirit  of  a  ” coarse  grained ”  examination  of  screw  pinches  in 
modern  pulselines  one  may  state  the  ” gas  bag'’  model  as  follows. 


The  differential  equaiton  for  the  relative  radial  motion  is 


a(a)  =  ( 


2l'£  iVgt  iivBt 


UB»  Ot  2  .  2  . 

--^rvs) 


with  ct  =  r(r)/r0,  0e  =  Te  +  (|)^ ,  and 

VS  =  \/{Z  +  1  )T{<dt)/mionrl  . 

The  axial  magnetic  frequencies  are  split  into  current  dependent,  compressing  exte¬ 
rior  (i)  and  initially  fixed,  retarding  interior  (ii)  components;  while  the  azimuthal 
magnetic  frequency  is  related  to  a  generalized  Bennett  current,  viz. 


ttNI(t) 

-  Cv/Y  ’  liUB‘  ~ 


4tt  mt0nn0rl 


vb*  = 


To  v7 irc2ml0nn0 


The  inductance  and  impedance  can  be  shown  to  be 


L  =  1-8  •  1021/p/ c2  (In  —  +  U  r2  -  r2)^]2)  [ni/]  , 

r  L 


Z  =  -  ( -  +  [2nN]2fr)  1.8  •  10 ulp/c2  (17]  , 

r 

Z0  =  3.257  •  10_6(/p/rrr2)(-^- )  — y  [0]  . 

1U  2  T  2 


Here  lp  is  the  pinch  length,  N  is  the  return  current  winding  number  (per  cm),  and 
rw  is  the  radius  of  the  ret  urn  current  cage.  The  factors  T  and  n0  are  the  pinch  mass 
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per  cm  and  mean  ion  density;  for  a  strict ly  self-similar  Gaussian  they  are  related 
by  T  =  motomn0iTa'10.  The  heating  equation  can  be  written 


-  -\i ^<e)  +  l02l~¥T^f 


+  z; 


I2Z 


po  2 

ct 2  77 


rad 


where  Zpo  =  Z0{  1.0375733  •  10  11  A/lp~f)  is  the  specific  impedance  (per  ion). 

Because  of  the  strong  radial  scaling  of  the  type  i  axial  magnetic  compression 
it  is  possible  for  this  model  to  fall  into  a  radiative  collapse  phase  toward  the  end  of 
an  implosion  cycle.  In  contrast  to  the  more  common  notion  of  refrigerative  collapse 
however,  the  compression  here  is  driven  from  the  exterior  axial  magnetic  field  and 
cannot  terminate  unless  either  the  current  is  depleted  from  the  pinch  due  to  high 
terminal  impedance  or  the  equation  of  state  changes  rapidly  to  stiffen  the  radial 
scaling  of  the  retarding  pressure  term.  Indeed  so  much  increase  in  density  (>1023 
ions  per  cm3 )  is  implied  for  this  class  of  solution  that  we  will  examine  its  relevance 
separately  and  here  seek  its  border  in  the  parameter  space  mainly  to  avoid  it. 


B.  A  Screw  Pinch  Load  in  the  Decade  Machine  Environment 

The  present  experiment  at  Maxwell  Laboratories  takes  it  roots  in  a  scheme  due 
to  Colgate6  for  stabilizing  a  pinch  without  having  to  compress  the  axial  field  which 
is  providing  the  stability.  By  twisting  the  return  current  cage  in  the  diode  front 
end  some  axial  magnetic  field  is  produced  by  the  full  driver  current  and  the  load 
can  he  stabilized.  This  means  to  a  fast  PRS  load  appears  to  suffer  coupling  losses 
if  driven  by  decade  class  machine  models. 

The  gas  bag  model  above  is  used  in  conjunction  with  the  transmission  line 
model  of  Ref.  1,  and  the  driver  is  chosen  as  a  model  of  the  Falcon  xlO  machine.  As 
the  windings  per  centimeter  are  increased  past  about  0.125  to  a  value  of  0.25,  the 
energy  delivery  to  the  model  load  falls  off  by  about  10%  and  will  continue  to  decay 
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because  of  the  added  inductance.  The  typical  implosion  scenario  of  such  a  screw 
pinch  load  in  concert  with  a  Falcon  type  inductive  store  is  shown  in  Figure  1(a). 
The  typical  performance  of  a  screw  pinch  load  as  a  function  of  the  N  value,  all  other 
parameters  being  equal,  is  shown  in  Figure  1(b).  The  use  of  a  high  winding  number 
for  stability  is  apparently  not  a  very  attractive  strategy  once  N  exceeds  0.125  or  so. 
The  act  ual  stability  requirement  is  not  known  for  this  particular  set  of  held  prohles, 
but  will  be  left  to  future  work. 
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Falcon  Driver  and  Gas  Bag  PRS  Screwpinch 


■ 


Figure  1(a) 


no 


b)  Ring  Plasma  Stagnation  Dynamics 

We  investigate  the  radiative  properties  and  dynamics  of  an  argon 
ring-puff  plasma  as  it  collides  with  and  stagnates  on  a  planar  iron 
surface.  Hydrogen  ring  plasmas  have  been  produced  experimentally 
with  velocities  exceeding  10s  cm/sec,  and  should  thus  be  capable  of 
generating  a  high  temperature,  high  density  plasma  on  impact  with  a 
target.  Such  a  plasma  would  efficiently  convert  the  kinetic  energy 
of  the  ring  plasma  into  multi-keV  radiation.  The  interaction  of  the 
ring  puff  with  the  target  is  simulated  with  a  ID  non-LTE 
multimaterial  radiation  hydrodynamics  model  which  self-consistently 
includes  the  effects  of  opacity  and  line  broadening.  This  study  has 
two  principal  goals:  to  determine  the  importance  of  wall  losses,  and 
to  characterize  the  radiation  produced  in  the  stagnation  process. 

The  puff  gas  was  modeled  in  ID  by  a  uniform  argon  plasma  moving 
at  uniform  velocity  u0  into  the  iron  wall  plasma.  This  is  a 
reasonable  approximation  for  timescales  less  than  D/u0 ,  where  D  is 
the  diameter  of  the  ring.  Argon  was  selected  as  the  gas-puff 
material,  since  hydrogen  (and  other  low-Z)  materials  become  fully 
stripped  in  a  high-velocity  stagnation,  and  cease  to  radiate 
efficiently.  By  modeling  the  interaction  in  this  way,  we  have 
emphasized  the  stagnation  and  radiation  physics,  and  ignored 
questions  relating  to  the  acceleration  of  the  ring-puff  and  its 
stability.  However,  we  view  these  latter  issues  as  technological, 
since  hydrogen  puffs  with  the  appropriate  characteristics  have 
already  been  produced. 
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A  benchmark  simulation  was  first  performed,  where  ionization, 
radiation  and  thermal  conduction  were  neglected.  This  idealized  case 
has  an  analytic  solution:  the  stagnation  results  in  the  formation  of 
a  layer  of  shocked  argon  (the  iron  remains  cold),  with  a  density 
ratio  of  4:1  for  large  Mach  numbers,  and  the  temperature  of  the 
shocked  argon  is  given  by  kT  =  1 /3 mj  u0  2  ,  where  is  the  mass  of  a 
gas-puff  ion,  for  the  case  of  no  ionization.  Thus,  for  the  idealized 
case  (no  wall  or  radiative  losses,  no  ionization,  etc.),  the 
stagnation  temperature  scales  as  the  square  of  the  puff  velocity.  A 
schematic  of  this  situation  is  given  in  Figure  1.  The  computer 
simulation  reproduced  the  analytic  model  fairly  accurately;  it  is  a 
severe  test  of  the  code,  since  the  gradients  are  very  large  and  the 
computational  mesh  must  be  severely  stretched  to  provide  adequate 
resolution  in  the  interface  region. 

Full  radiation-hydrodynamic  simulations  with  detailed  atomic 
models  for  iron  and  argon  were  performed  for  several  puff  velocities: 
u0=lxl08,  u0 =2x10®  and  u0 =3x10® cm/sec.  In  each  of  these  cases,  a 
hot  argon  layer  is  created  initially  at  the  interface  with  the  wall 
plasma,  which  grows  in  size  with  time.  A  weak  shock  wave  propagates 
into  the  wall  plasma,  and  energy  transport  by  radiation  and  thermal 
conduction  heats  the  iron  plasma  and  produces  an  ablation  into  the 
argon  puff.  An  energy  history  for  the  2x10® cm/sec  case  is  shown  in 
Figure  2.  Initially,  most  of  the  energy  in  the  system  is  in  the 
kinetic  energy  of  the  argon  puff.  As  the  stagnation  onto  the  wall 
proceeds,  kinetic  energy  is  converted  into  thermal  energy,  and  the 
heated  plasma  ionizes  (ionization  energy)  and  radiates  (radiated 
energy).  By  about  0.45  nsec,  most  of  the  kinetic  energy  is 
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thermalized,  and  the  hot  plasma  begins  to  expand.  In  the  expansion 
phase,  thermal  energy  is  converted  into  kinetic,  and  ionization 
energy  is  converted  into  thermal  (and  radiation)  through 
recombination.  The  radiated  energy  is  lost  to  the  system. 

In  Figure  3.  the  hydrodynamic  variables  (density  [g/cm3], 
temperature  [eV],  pressure  [ergs/cm3 xlO' 1 2  ] ,  velocity  [cm/nsec  xlO'2] 
and  radiative  cooling  [ergs/cm3 -nsec  xlCr14])  are  plotted  versus 
displacement  in  centimeters  for  the  2xl08 cm/sec  case  at  0.33  nsec, 
just  before  the  minimum  in  kinetic  energy.  It  car.  be  seen  that  the 
frontside  of  the  iron  plasma  is  hot  (about  2.5  keV),  and  blowing  off 
into  the  argon  plasma.  The  radiative  cooling  peaks  just  behind  this 
blowoff  plasma,  where  the  density  and  temperature  are  both  moderately 
high.  The  temperature  in  the  argon  puff  is  nearly  uniform  (about  3-0 
keV)  at  this  time  because  of  thermal  conduction  and  radiation 
transport.  The  emission  spectrum  at  this  time  is  shown  in  Figure  4. 
This  spectrum  and  the  radiation  data  in  the  tables  described  below 
were  post-processed  using  extended  atomic  models  (more  detailed  than 
required  for  energetics).  The  spectrum  is  characteristic  of  a  high 
density  plasma;  continuum  radiation  dominates  the  line  emission.  In 
contrast,  wire  plasmas  and  ordinary  puff-gas  plasmas  generally  reach 
lower  densities  than  are  achieved  in  the  present  configuration,  and 
are  less  dominated  by  continuum  emission. 

For  puff-gas  velocities  substantially  smaller  than  2x10® cm/sec, 
much  of  the  energy  radiated  above  2  keV  comes  from  argon  K-shell 
lines  (principally  from  the  transitions  n=2  to  n=l  of  Ar  XVIII  and 
ls2p{lP)  to  Is2  of  Ar  XVII);  however,  in  the  present  case,  most  of 


the  multi-keV  energy  is  produced  by  bound-free  processes  in  the  hot, 


dense  iron  plasma  (principally  recombination  into  lithium-like 
states).  A  breakdown  of  the  instantaneous  radiated  powers  at  0.33 
nsec  by  element  and  by  type  (line,  recombination  and  bremsstrahlung) 
is  given  in  Table  I.  It  can  be  seen  that  the  total  radiated  power 
from  the  iron  plasma  is  almost  an  order  of  magnitude  higher  than  that 
from  the  argon;  this  ratio  increases  with  puff  velocity.  The 
radiation  is  broken  down  into  several  energy  ranges  in  Table  II.  At 
0.33  nsec,  the  power  peaks  in  the  2-4  keV  range;  however,  there  is 
comparable  power  radiated  in  the  1-2  keV,  2-4  keV  and  >4.0  keV  bands. 
The  largest  contribution  to  the  power  above  1  keV  is  bound-free 
radiation  from  the  relatively  dense  iron  plasma,  as  described  above. 

A  temperature  history  plot  (log10  temperature  contours  versus 
position  and  time)  is  given  in  Figure  5  Tor  the  2x10® cm/sec  case. 
Each  horizontal  slice  represents  a  .^tc-efic  time,  and  the  bottom  edge 
corresponds  to  the  init  al  conditions  (7-5eV  argon  and  0.7eV  iron). 
The  void  at  the  rig!  t  side  of  the  plot  corresponds  to  the  outer  edge 
of  the  ring-puff,  which  moves  inwards  until  the  stagnation  is 
complete  at  about  0.45  nsec,  and  then  moves  outwards  in  the 
expansion  phase.  The  peak  temperature  calculated  in  the  simulation 
is  a  little  more  than  3  keV,  whereas  the  stagnation  temperature  for 
an  ideal  gas  (kTs  =  1 /3 m1 u0 2 )  is  533*3  keV.  Hovever,  the  argon 
ionizes  to  the  H  and  He-like  state  ( 1 6- 17  times  ionized),  and  much  of 
the  iron  ionizes  to  the  Li  and  Be-like  state  (22-23  times  ionized) . 
Therefore,  most  of  the  kinetic  energy  goes  into  heating  the 
electrons.  Additionally,  radiation  and  thermal  conduction  carry 
energy  into  the  dense  wall  plasma  and  out  into  the  unshocked  puff- 
gas.  Finally,  energy  is  expended  in  ionizing  the  plasma  and  energy 
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is  radiated  away. 


By  about  0.1  nsec  the  puff  gas  has  become  nearly  isothermal 
(about  3  keV).  The  stagnation  phase  ends  at  about  0.45  nsec;  at  this 
time,  the  argon  is  beginning  to  cool.  The  front  side  of  the  iron 
becomes  warm  and  begins  to  expand  into  the  gas  puff  by  0.1  nsec.  The 
iron/argon  interface  is  shown  in  Figure  6.  The  heating  and  expansion 
continue  throughout  the  stagnation;  when  the  gas  puff  begins  to 
expand  and  cool,  the  front  side  of  the  target  also  begins  to  cool. 
However,  energy  continues  to  be  transported  into  the  iron  plasma,  and 
this  energy  is  radiated  away  at  a  rate  that  maintains  the  rearside 
temperature  near  100  eV  from  about  0.3  to  beyond  1.2  nsec. 

A  history  plot  of  the  radiated  power  is  given  in  Figure  7- 
There  are  two  distinct  radiating  regions:  one  at  the  frontside  of  the 
iron  plasma,  and  the  other  in  the  argon  plasma  just  ahead  of  the 
iron/argon  interface.  Net  radiation  loss  is  plotted;  there  are 
regions  that  are  net  absorbers  of  radiation.  The  radiation  cooling 
is  the  greatest  in  places  where  temperature  and  density  are 
simultaneously  large.  The  cooling  peaks  at  about  the  same  time  that 
the  kinetic  energy  is  a  minimum,  and  it  declines  (maintaining  two 
spatially  distinct  regions)  as  the  plasma  expands  and  cools. 

Some  idea  of  the  complexity  of  the  puff-target  stagnation  can  be 
seen  from  the  history  plot  of  mass  density  for  the  2xl08cm/sec  case, 
shown  in  Figure  8.  Substantial  puff-gas  compression  occurs  initially 
near  the  iron  interface,  but  expansion  takes  place  just  behind  this 
region  as  a  consequence  of  local  heating  (due  to  radiation  transport 
and  thermal  conduction).  Eventually,  the  outer  edge  of  the  argon 
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plasma  becomes  rarefied  because  of  energy  transport  deep  into  the 
gas-puff.  The  ablation  of  the  iron  target  can  also  be  discerned,  as 
well  as  alternate  compression  and  rarefaction  near  the  interface  as 
localized  "bouncing"  occurs .  The  density  in  the  bulk  of  the  iron 
plasma  monotonically  decreases  due  to  heating  and  subsequent  ablation 
into  the  gas  puff. 

Temperature  history  plots  are  given  in  Figures  9  and  10  for 
the  lxl08cm/sec  and  3*10® cm/sec  cases,  respectively.  The  most 
obvious  trend  which  can  be  seen  from  the  temperature  histories  is  the 
dependence  of  the  "bounce"  time  on  u0  the  initial  gas-puff  velocity. 
Another  obvious  trend  is  the  increase  of  temperature,  both  in  the 
puff  and  in  the  iron  plasma,  with  increasing  u0 .  Peak  temperatures 
reach  about  1.8,  3-0  and  4.0  kev  for  the  1.0,  2.0  and  3 .0x10® cm/sec 
cases.  Finally,  it  can  be  seen  that  transport  of  energy  into  the 
puff  becomes  more  efficient  at  higher  puff  velocities.  A  nearly 
isothermal  puff-gas  is  produced  almost  immediately  at  3*10® cm/sec;  it 
requires  about  a  half  a  nanosecond  for  thermal  gradients  in  the  puff 
to  be  correspondingly  reduced  at  1x10® cm/sec. 

It  has  been  confirmed  that  the  stagnation  of  a  high  velocity 
ring-puff  plasma  on  a  planar  solid  density  target  can  generate  a  high 
density,  high  temperature  plasma,  which  can  efficiently  generate 
multi-keV  radiation.  However,  several  questions  remain.  These 
include  the  importance  of  2D  effects  which  were  ignored  in  this 
study,  and  the  importance  of  non-fluid  behavior  at  the  puff-target 
interface  (high-velocity  ring  plasma  ions  may  penetrate  the  target 
and,  in  extreme  cases,  the  interaction  may  resemble  a  deposition 
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rather  than  a  stagnation  if  the  B- field  insulation  breaks  down) .  In 
addition,  the  adequacy  of  the  atomic  models  employed  in  this  study 
should  be  further  investigated.  In  particular,  the  adequacy  of  the 
structure  in  the  lower  ionized  states  of  iron  should  be  scrutinized. 
Finally,  there  may  be  magnetic  {or  electric)  fields  produced  in  the 
acceleration  of  the  ring  puff,  or  in  its  interaction  with  the  target, 
which  could  affect  the  dynamics.  These  complications  are  currently 
under  investigation. 
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TABLE  I 

Instantaneous  radiated  powers  (W/cm2 ) 
[2x10s  cm/sec  case  at  0.33  nsec] 


Iron 

Argon 

Total 

Bound-bound 

1. 52x10* 3 

2.34xl012 

1.75x10* 3 

Bound- free 

5.22x10* 3 

5.49x10* 2 

5.77xl013 

Free- free 

4. 38X101 2 

1.07xl012 

5.45x10* 2 

Total 

7.18x10* 3 

8.90xl012 

8.07xl013 

TABLE  II 

Instantaneous 

radiated  powers 

(W/cm2 ) 

[2xl08 cm/sec  case  at  0.33 

nsec] 

Bound-Bound 

Bound -Free 

Free-Free 

Total 

0.0  -  0.1  keV 

2.73x109 

3.92x10® 

5.27x10** 

5.30X1011 

0.1  -  0.2  keV 

9.02xl018 

3.30X109 

3.85x10** 

4.79xl031 

0.2  -  0.4  keV 

5-57x10** 

5.03xl011 

5-58x10** 

1.62xl012 

0.4  -  1.0  keV 

8.82X1011 

6.39xl012 

1.19xl012 

8.47x1012 

1.0  -  2.0  keV 

l.OlxlO13 

1.04xl013 

1.21xl012 

2.17xl013 

2.0  -  4.0  keV 

2.13xl012 

2.80x10* 3 

1.06xl012 

3.12xl013 

above  4 . 0  keV 

3-77x10* 2 

1.25xl013 

5.16x10** 

1.67xl013 

Total 

1.75x10* 3 

5.77xl013 

5.45x1012 

8.07xl013 
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LOG  TEMPERATURE  VS  TIME 
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c)  Evolution  of  An  X-Pinch  Plasma 
I.  Introduction 

The  hot  dense  plasma  that  results  from  the  passage  of  large  currents 

through  a  Z-pinch  plasma  has  attracted  considerable  attention  in  recent 

years  as  a  radiation  source  for  x-ray  laser  research,  x-ray  lithography 

applications,  broadband  x-ray  flashlamp  studies,  and  spectroscopy 

experiments^.  The  majority  of  experiments  and  theoretical  investigations 

have  been  performed  for  the  standard  or  conventional  load  designs  such  as 

single  and  multiple  wires,  gas  puffs,  metal  foils,  and  more  recently,  metal 

vapors.  With  few  exceptions  almost  all  these  loads  exhibit  unstable 

behavior  during  the  implosion  on  axis.  By  altering  the  configuration  to 

include  an  inner  axial  magnetic  field  or  an  inner  cylindrical  plasma 

localized  along  the  axis,  the  m=0  instability  can  be  inhibited.  Both  of 

these  schemes  tend  to  stabilize  the  implosion  and  have  the  effect  of  a 

stagnation  post.  However,  the  K-shell  radiative  efficiency  of  such 

configurations  is  somewhat  reduced  from,  say,  a  single  wire  load  even  though 

the  radiative  cooling  rate  summed  over  all  atomic  shells  is  virtually  the 

2 

same  for  identical  generator  and  load  conditions  .  For  example,  a  single 
aluminum  wire  upon  heating,  first  expands  and  subsequently  runs-in  and 
assembles  on  axis  giving  rise  to  a  series  of  bright  or  hot  spots.  These  hot 
spots  are  distributed  along  the  axis  and  radiate  intense  short  bursts  of  x- 
rays  corresponding  to  temperatures  and  densities  which  are  higher  than  the 
corresponding  values  associated  with  the  bulk  plasma.  The  assembled  plasma 
achieves  greater  compressions  leading  to  higher  densities  and  smaller  radii 
in  the  pinched  regions  than  do  the  other  types  of  load  configurations  such 
as  wire  arrays  or  gas  puffs.  Since  the  radiative  cooling  rate  varies 
roughly  as  the  density  squared  (for  free  streaming  optically  thin  plasmas) 
multiplied  by  the  characteristic  radius  of  the  spot  squared,  the  K-shell 
radiative  efficiency  is  also  greater  for  single  wires,  all  else  being  equal. 

It  is  this  behavior  that  has  attracted  our  attention  to  the  x-pinch 
configuration.  Since  the  m=0  sausage  instability  occurs  as  a  natural 
consequence  of  the  implosion  dynamics  for  conventional  loads  producing 


enhanced  temperatures  and  densities,  a  logical  extension  in  load  design 
would  be  to  incorporate  these  features  initially  in  order  to  take  advantage 
of  the  higher  temperatures  and  densities  and,  most  importantly,  improved 
radiative  efficiency.  In  fact,  the  x-pinch  configuration  has  been 
investigated  experimentally  for  currents  at  the  100  kiloampere  level  and  was 
found  to  behave  as  anticipated;  the  pinched  region  or  neck  was  clearly 
distinguished  in  pinhole  photographs  as  a  localized  spot  of  bright  emission 
of  thermal  origin  for  the  most  part  .  The  measured  size  of  the  spot,  based 
on  small  pinhole  estimates,  was  about  10  pm.  The  spectrum  revealed  K-shell 
emissions  from  aluminum  and  L-shell  emissions  from  nickel.  Multiple  copper 

x-pinch  loads  at  higher  current  levels  has  also  demonstrated  K-shell 

.  .  4 
emission  . 

In  this  paper  we  investigate  the  dynamics  and  radiative  properties  of 
an  array  of  thin  copper  wires  which  are  twisted  in  the  center  and  configured 
to  form  an  x-pinch.  The  dynamics  are  described  by  a  2-D  non-LTE  radiation 
magneto-hydrodynamics  model  (PRISM)  coupled  to  an  external  driving  current. 
The  symmetry  imposed  by  the  r-z  coordinates  (i.e.,  the  0  direction  is 
ignored)  requires  that  we  model  the  problem  after  it  has  been  assumed  that 
the  plasmas  from  the  wires  have  merged  and  have  the  form  rather  like  that  of 
a  hollow  cone.  Individual  wires  in  the  array,  therefore,  cannot  be  modeled 
as  separate  entities.  The  focus  of  this  investigation  is  on  the  dynamic 
evolution  and  morphology  of  the  plasma  rather  than  the  determination  of  an 
optimal  design  configuration  for  enhanced  radiation  efficiency  or  x-ray 
laser  scenarios.  It  is  our  intent  to  gain  a  better  understanding  of  these 
plasmas  in  order  to  compare  and  contrast  the  results  in  the  future  with 
experimental  observations. 

II.  Model 

The  basic  hydrodynamic  quantities  of  mass,  momentum,  and  total  energy 
are  transported  using  a  Lagrangian  scheme  on  a  two-dimensional  mesh.  A 
remap  phase  using  parabolic  interpolation  is  then  employed  to  update  the 
variables  on  the  original  mesh.  The  dynamic  evolution  of  the  system  is 
represented  as  solutions  to  the  following  equations: 
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To  these  equations  must  be  added  Maxwell's  equations,  viz. 

V  x  B  .  —  J 
~  ~  c 


and  Ohm's  law, 


.  9B 

V  x  E  =  -  - 

c  at 


E=qj--uxB 
-  ~  c  ~ 


where  P.  =  ion  (electron)  pressure,  C  .  is  the  electron-ion  energy  exchange 
term,  e.  =  ion  (electron)  specific  energy,  p  =  density,  V  =  1/p,  u  =  (u  , 
uz)  =  velocity,  q  =  heat  flux,  [q  =  -  K^V  T  -  K^V  T  -  KA  (eB  x  V  T)  ], 

q  =  ohmic  heating,  and  P  .  radiation  power.  The  resistivity  model  of 

5  raa  6 

Spitzer  is  used  along  with  the  relevant  Braginskii  transport  coefficients  . 

The  total  current,  I,  is  specified  as  a  function  of  time  and  sets  the  outside 

boundary  condition  on  the  magnetic  field  according  to  B0  =  21/rc.  All  other 

V 

symbols  have  their  conventional  meaning.  The  pressure  is  determined  from  the 
equation  of  state,  viz. 
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where 


=  P<p’eint> 


0  =0  +0 
int  thermal  ionization 

For  example,  in  the  case  of  a  perfect  gas 

P  *  p  thermal  ' 

In  our  case  (p,  )  are  obtained  from  the  above  equations  and  are  then  used  to 

evaluate  P  and  T. 

The  atomic  level  populations  in  the  plasma  are  determined  by  a  set  of  rate 
equations  of  the  form: 


Z 

j 


V.  .f . 
Ji  3 


£  V..f., 
j  1J  1 


(1) 


where  f^  is  the  fractional  population  of  atomic  level  i,  and  is  the  net 
reaction  rate  describing  the  transition  from  initial  state  j  to  final  state  i. 
An  equation  of  this  type  is  constructed  for  each  of  the  atomic  levels  included 
in  the  model.  For  sufficiently  dense  plasmas,  the  effective  populating  and 
depopulating  rates  are  generally  fast  compared  with  the  hydrodynamic  response. 
Under  these  circumstances,  an  equilibrium  assumption,  which  involves  dropping 
the  explicit  time  dependence  in  Eq.  (1),  can  be  justified.  The  plasma  is  then 
said  to  be  in  collisional-radiative  equilibrium  (CRE),  allowing  the  plasma 
ionization  state  to  respond  instantaneously  to  changes  in  hydrodynamic 
quantities . 


The  rate  coefficients  that  are  used  to  calculate  the  populating  and 
depopulating  rates,  U^,  are  calculated  using  various  atomic  calculational 
methods.  The  processes  included  in  this  calculation  and  the  methods  used  in 
calculating  the  corresponding  rate  coefficients  are  summarized  elsewhere^.  Once 
the  set  of  rate  equations  (including  the  optical  pumping  from  the  radiation 
field)  has  been  solved  for  the  level  populations,  f^,  the  electron  density  can 
be  calculated. 
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Ne  =  Z.  Z.  £•  Nr 


where  Z..  is  the  ionic  charge  of  level  i  and  is  the  total  ion  density. 

The  ionization  and  excitation  energy  can  also  be  calculated  by 

eT  =  Z  X.  f.  NT, 

Ionization  ^  1  1  I 

where  X.  is  the  energy  of  level  i,  measured  from  the  ground  state  of  the  neutral 
atom.  The  copper  model  contains  all  the  ground  states  from  the  neutral  to  the 
fully  ionized  state.  In  addition,  there  is  excited  state  structure  in  the  H-, 
He-,  Li-,  and  Be-like  ionization  stages,  containing  10  emission  lines  selected 
for  their  energetic  and  diagnostic  value. 

The  local  rate  of  energy  change  in  zone  j,  a  result  of  radiative  processes, 
is  given  by 

■  -  f  FPi 

where  is  the  rate  of  energy  loss  in  zone  j  resulting  from  a  discrete 

radiative  process  (or  frequency  group)  P,  for  each  bound-bound,  bound-free,  and 
free-free  process. 


III.  Results  and  Discussion 


In  the  numerical  simulations,  we  have  assumed  that  the  wires  are  vaporized 
and  a  plasma  has  been  created  which  is  symmetric  around  the  z-axis  with  a 
conical  cross-section.  With  this  simplification,  the  2-D  simulation  is  carried 
out  in  the  r-z  plane.  We  have  chosen,  somewhat  arbitrarily,  a  midplane  density 
which  is  bounded  at  0.075  cm  in  both  r  and  z.  The  cone  thickness  is  tapered  to 
a  thickness  of  0.02  cm  at  r=0.3  cm  and  z=1.0  cm.  The  same  conditions  exist  in 
the  other  three  quarter  planes.  The  total  length  in  the  z-direction  (i.e.,  the 
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diode  length)  is  then  2  cm.  We  have  examined  cases  in  which  either  the  mass  per 
unit  length  along  z  or  the  density  within  the  cone  is  taken  to  be  uniform.  The 
latter  case  can  be  used  to  model,  for  instance,  an  etched  or  tapered  wire  array. 
In  this  section  we  consider  only  cases  with  an  initial  constant  density.  The 
density  and  radius  are  chosen  to  correspond  to  an  initial  array  of  copper  wires 
twisted  at  the  center.  Outside  of  this  region,  the  density  is  set  to  a  floor 
which  is  typically  less  than  12  of  the  central  interior  value. 

A.  Hydrodynamic  Simulations 

We  first  examine  the  phenomenology  of  an  x-pinch  implosion  with  a  simulation 

in  which  there  are  no  radiative  loss  terms.  Figure  1  shows  the  initial  density 

contours  over  the  half-plane  z=0  to  z=1.0  cm.  In  this  example  we  have  chosen  an 

-2  3 

initial  uniform  density  of  4  x  10  gm/cm  and  an  initial  electron  temperature 
of  30  eV.  The  initial  ion  internal  energy  is  chosen  to  be  a  factor  of  ten  below 
the  electron  internal  energy.  An  external  current  is  applied  as  a  boundary 

condition  on  rBQ  and  is  chosen  to  rise  to  5  MA  in  25  nsec.  In  this  simulation 

“  8 
we  have  used  the  equation  of  state  tables,  SESAME,  for  copper  .  The  mesh 

resolution  is  124  in  the  z-direction  and  74  in  the  r-direction  on  the  quarter- 

plane.  The  values  on  the  opposite  quarter  plane  in  the  contour  plots  are 

determined  from  symmetry  arguments. 

Figure  2  shows  the  density  contours  at  10,  14,  18,  and  26  nsec.  The 

calculation  is  terminated  after  26  nsec  when  plasma  has  been  pushed  up  against 
the  top  wall  which  is  located  at  1  cm.  At  10  nsec,  the  density  contours 
indicate  that  there  is  plasma  motion  up  the  z-axis  as  a  result  of  the  stagnation 
of  the  inside  part  of  the  cone  on  axis.  By  14  nsec  the  peak  density  has  been 
localized  in  a  spot  which  is  above  the  r-axis.  This  spot  has  moved  further  up 
the  axis  by  18  nsec,  and  further  still  at  the  last  frame,  26  nsec.  Figure  3 
displays  logarithmic  electron  temperature  and  density  curves  along  the  z-axis. 
The  frames  are  at  the  same  times  as  shown  in  Fig.  2.  The  peak  density  is  seen 
to  move  up  the  z-axis.  A  cavity  or  low  density  region  evolves  behind  the  peak 
density  due  to  material  being  pushed  away  from  the  x-point.  The  back  edge  of 
this  spot  (as  seen  in  Fig.  2(c))  has  reached  an  electron  temperature  of  1  keV 
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while  the  bulk  of  the  spot  has  a  temperature  of  approximately  100  eV.  This 
gradient  eventually  diffuses  through  the  high  density  region  and  the  entire  spot 
is  roughly  200  eV  at  18  nsec.  In  the  last  frame  of  this  Figure,  the  temperature 
has  again  risen  above  the  keV  level  in  the  underdense  cavity  behind  the  peak 
density.  This  is  a  result  of  material  being  compressed  into  this  region  and  the 
plasma  attempting  to  establish  equilibrium.  A  low  density  foot  has  arrived  at 
the  upper  wall  at  this  time  and  the  simulation  is  subsequently  terminated. 

B.  Radiation  Hydrodynamic  Simulations 

We  have  investigated  different  mass  loadings  or  equivalently,  wire  diameters 
for  the  generic  design  discussed  above  (i.e.,  a  tapered  cone  which  intersects 
the  r-  and  z-axis  at  0.075  cm,  and  a  0.02  cm  width  at  r  =  0.3  cm  and  z  =  1.0  cm) 
but  with  radiation  loss  terms  included.  In  these  simulations,  we  make  use  of  an 
atomic  model  for  Copper  which  contains  all  the  ground  states  as  well  as 
structure  for  the  H-,  He-,  Li-,  and  Be-like  ionization  stages.  Ten  transitions 
(lines)  within  each  of  these  stages,  as  well  as  the  bound  free  edges  for  all 
excited  and  ground  states  are  carried.  The  radiated  power  which  includes  the 
bound-bound,  bound-free,  and  free-free  components  is  then  used  to  calculate  a 
total  radiated  power  for  each  computational  zone.  Internal  energy  is 
partitioned  into  ionization  and  thermal  components  so  that  electron  temperature 
and  pressure  are  self-consistently  calculated.  These  quantities  are  then  used 
instead  of  those  obtained  from  the  SESAME  tables  as  in  the  previous  section. 
The  ions  are  treated  as  a  fluid  with  y=5/3.  The  large  amount  of  memory  and 
computational  time  required  for  this  model  dictates  that  we  limit  ourselves  to  a 
mesh  of  size  35  x  60.  Here,  we  do  not  consider  radiation  transport  with  the 
proviso  that  the  maximum  power  output  from  a  given  zone  is  limited  to  a  fraction 
of  the  Blackbody  limit.  This  limiter  is  usually  set  to  0.1. 

_2 

Ve  have  considered  four  different  mass  densities  in  the  cone:  2  x  10  ,  4  x 
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10  ,  8  x  10  ,  and  12  x  10  gm/cm  .  In  the  following  discussion  we  again 

-2  3 

concentrate  on  the  case  where  the  initial  density  is  4  x  10  gm/cm  .  Figure  4 
shows  density  contours  at  8,  14,  16,  and  18  nsec.  These  plots  show 

qualitatively  the  same  behavior  as  seen  in  Fig.  2  with  motion  of  a  dense  spot 
along  the  z-axis.  However,  multiple  spot  formation  can  be  seen  near  the  origin 


as  evidenced  by  the  contours  at  16  nsec.  In  the  SESAME  hydro  only  simulation, 
the  dense  spot  moved  up  the  axis  and  left  a  cavity  behind.  In  this  case,  a 
dense  region  reforms  near  the  origin  so  that  by  18  nsec  there  are  two  distinct 
dense  regions  which  can  be  seen.  Other,  less  dense  regions  have  also  formed 
along  the  outside  of  the  cone.  The  foot  in  front  of  the  dense  region  also 

extends  further  up  the  z-axis  than  in  the  hydro  only  run  and  contacts  the  upper 
wall  shortly  after  18  nsec. 

The  logarithm  of  the  radiated  power  is  shown  in  Fig.  5  at  8,  14,  16,  and  18 
nsec.  In  the  first  frame,  8  nsec,  the  highest  levels  of  radiated  power  occurs 
near  the  outside  of  the  pinch  plasma  in  the  region  of  the  x-point.  Radiation  is 
also  coming  from  material  which  is  being  pushed  up  the  axis.  The  low  outside 
material  density  combined  with  the  high  current  density  in  this  region  tends  to 
produce  ohmic  heating  levels  which  initially  dominate  the  radiation.  Although 
current  is  continuously  diffusing  towards  the  axis,  the  current  density  remains 
highest  on  the  outside  while  the  total  current  is  increasing.  This  contributes 
substantial  ohmic  heating  and  radiative  power.  As  the  pinch  develops  further, 
at  14  and  16  nsec,  the  radiated  power  is  maximized  on  the  outside  of  the  dense 
spots  seen  in  the  density  contours  of  Fig.  4.  At  18  nsec,  the  upper  and  lower 
dense  regions  of  Fig.  4(d)  are  seen  to  be  radiating  most  strongly.  The  location 
of  the  upper  region  is  at  approximately  the  same  position  as  the  upper  dense 

region  shown  in  Fig.  2(d)  and  Fig.  4(d).  The  bottom  region  represents  a  large, 

radiating  hot  spot  located  at  the  x-point. 

The  total  radiated  line  energy  as  a  function  of  initial  x-pinch  density  is 
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shown  in  Fig.  6.  At  the  two  lowest  densities,  2  x  10  gm/cm  and  4  x  10 

3 

gm/cm  (which  would  correspond  to  an  array  of  four,  41.5  ym  and  83  urn  diameter 

copper  wires  respectively)  the  yield  is  roughly  about  a  kilojoule  in  all  lines 

2  3 

and  falls  to  about  150  J  at  the  highest  density  of  12  x  10  gm/cm  .  The  yield 

for  the  lines  above  8  keV  in  energy  (i.e.,  K-shell  lines)  are  45  J  for  the  2  x 

-23  -23 

10  gm/cm  and  13  J  for  the  4  x  10  gm/cm  case.  This  drops  below  10  J  for 

the  higher  densities  or  equivalently,  larger  wire  diameter  arrays.  In  contrast 

to  this,  the  bound-free  continuum  contribution  (i.e.,  from  all  ground  states  as 

_o 

well  as  the  excited  states)  to  the  radiative  yield  was  32  kJ  for  the  4  x  10 

3 

gm/cm  simulation.  All  yield  values  were  calculated  up  to  18  nsec  when  the 


137 


calculations  were  terminated  because  of  material  interacting  with  the  upper 

wall.  It  should  also  be  recognized  that  the  yield  values  given  above  and  shown 

in  Fig.  6  are  approximate  since  initial  starting  conditions,  x-pinch 

configuration,  mass  loadings,  etc.  can  change  these  values.  For  instance, 

_2 

lowering  the  initial  ion  internal  energy  by  a  factor  of  ten  in  the  2  x  10 

3 

gm/cm  case  produces  a  slightly  harder  stagnation  and  increased  the  K-shell 
yield  by  about  50%  although  the  total  yield  decreased  somewhat. 

-2  3 

The  radiation  power  spectrum  from  the  4  x  10  gm/cm  x-pinch  simulation  at 
18  nsec  is  shown  in  Fig.  7.  The  energy  scale  runs  from  200  eV  to  12  keV.  The 
logarithm  of  the  power  runs  over  six  orders  of  magnitude.  Line  broadening  due 
to  ion  thermal  motion  has  been  included  in  this  Figure.  However,  no 
instrumental  broadening  has  been  assumed.  The  dominant  lines  below  2  keV  are 
from  the  Li-  and  Be-like  stages.  The  weaker  higher  energy  lines  from  the  K- 
shell  (He-a,  Ly-a,  He-g,  and  Ly-g)  are  located  around  10  keV. 

IV.  Conclusion 

The  evolution  of  an  x-pinch  has  been  studied  from  the  context  of  a  single 
generic  configuration.  Ve  concentrated  on  the  dynamics  of  the  x-pinch  and  not 
on  an  optimization  study.  Ve  have  found  that  an  array  of  wires  twisted  to  form 
an  x-pinch  results  in  the  formation  of  a  hot  spot  or  jet  which  moves  up  the  z- 
axis.  In  the  case  which  was  presented  here,  a  hot  spot  was  also  produced  near 
the  x-point.  Radiated  power  was  concentrated  in  the  moving  spot  as  well  as  in 
the  central  hot  spot.  Some  Cu  K-shell  emission  (<100  J)  was  also  produced. 

A  more  exhaustive  study  will  have  to  be  conducted  in  order  to  determine  the 
importance  of  changing  the  current  pulse,  wire  loading  (e.g.,  tapering),  wire 
mass,  wire  material,  and  x-point  angle  to  the  improvement  or  degradation  of  the 
results  reported  on  here. 
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Figure  Captions 
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1.  Initial  configuration  for  an  initially  uniform  4  x  10  gm/cm  x-pinch. 

2.  Density  contours  at  various  times.  Eight  equally  and  linearly  spaced 
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contours  are  shown,  (a)  10  nsec  and  peak  density  of  8.7  x  10  gm/cm  . 

-1  3 

(b)  14  nsec  and  peak  density  of  3.45  x  10  gm/cm  .  (c)  18  nsec  and  peak 

-13  -1 

density  of  1.85  x  10  gm/cm  .  (d)  26  nsec  and  peak  density  of  3  x  10 

.3 

gm/cm  . 

3.  Logarithmic  temperature  and  density  curves  along  the  z-axis.  (a)  0 
nsec,  (b)  10  nsec,  (c)  14  nsec,  (d)  18  nsec,  and  (e)  26  nsec. 

_2 

4.  Density  contours  as  in  Fig.  2  but  for  an  initially  uniform  4  x  10 

3 

gm/cm  x-pinch  simulation  in  which  radiation  loss  terms  are  included,  (a) 
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8  nsec  and  peak  density  of  6.25  x  10  gm/cm  .  (b)  14  nsec  and  peak 
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density  of  7.4  x  10  gm/cm  .  (c)  16  nsec  and  peak  density  of  1.7  x  10 

3  -2  3 

gm/cm  .  (d)  18  nsec  and  peak  density  of  8  x  10  gm/cm  . 

3 

5.  Logarithmic  radiation  power  (ergs/sec-cm  )  contours  for  the  times  shown 

in  Fig.  4.  Each  plot  ranges  over  8  contours  with  3  contours  per  decade. 

21  3  22 

The  highest  possible  value  in  (a)  is  4.6  x  10  ergs/sec-cm  and  4.6  x  10 

3 

ergs/sec-cm  in  (b),  (c),  and  (d). 

6.  Total  line  radiation  as  a  function  of  initial  density. 

7.  Instantaneous  radiation  power  (ergs/sec-keV)  spectrum  as  a  function  of 
photon  energy  at  18  nsec. 
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(a) 

Fig.  2  —  Density  contours  at  various  times.  Eight  equally  and  linearly  spaced  contours  are  shown, 
(a)  10  nsec  and  peak  density  of  8.7  x  10“ 2  gm/cm3.  (b)  14  nsec  and  peak  density  of 
3.45  x  10-1  gm/cm3.  (c)  18  nsec  and  peak  density  of  1.85  x  10" 1  gm/cm3  (d)  26  nsec  and  peak 
density  of  3  x  10" 1  gm/cm3. 
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Fig.  2  (Continued)  —  Density  contours  at  various  times.  Eight  equally  and  linearly  spaced  contours 
are  shown,  (a)  10  nsec  and  peak  density  of  8.7  x  10" 2  gm/cm3.  (b)  14  nsec  and  peak  density  of 
3.45  x  10" 1  gm/cm3.  (c)  18  nsec  and  peak  density  of  1.85  x  10" 1  gm/cm3  (d)  26  nsec  and  peak 
density  of  3  x  10  1  gm/cm3. 


Fig.  2  (Continued)  —  Density  contours  at  various  times.  Eight  equally  and  linearly  spaced  contours 
are  shown,  (a)  10  nsec  and  peak  density  of  8.7  x  10“ 2  gm/cm5.  (b)  14  nsec  and  peak  density  of 
3.45  x  10” 1  gm/cm5.  (c)  18  nsec  and  peak  density  of  1.85  X  10“ 1  gm/cm5  (d)  26  nsec  and  peak 
density  of  3  x  10“ 1  gm/cm5. 


44 


:  10°  p  (gm/crrr),  Te  (eV) 


(a) 

Fig.  3  —  Logarithmic  temperature  and  density  curves  along  the  z-axis. 
(a)  0  nsec,  (b)  10  nsec,  (c)  14  nsec,  (d)  18  nsec,  and  (e)  26  nsec. 
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Fig.  3  (Continued)  —  Logarithmic  temperatun 
(a)  0  nsec,  (b)  10  nsec,  (c)  14  nsec, 
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(d) 


Fig.  3  (Continued)  —  Logarithmic  temperature  and  density  curves  along  the  z-axis. 
(a)  0  nsec,  (b)  10  nsec,  (c)  14  nsec,  (d)  18  nsec,  and  (e)  26  nsec. 
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Fig.  3  (Continued)  —  Logarithmic  temperature  at 
(a)  0  nsec,  (b)  10  nsec,  (c)  14  nsec,  (d) 
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simulation  in  which  radiation  loss  terms  are  included.  8  nsec  and  peak  density  of 


6.25  x  10  2  gm/cm5.  (b)  14  nsec  and  peak  density  of  7.4  x  10  2  gm/cms.  (c)  16  nsec  and  peak 
density  of  1.7  x  10-1  gm/cm5.  (d)  18  nsec  and  peak  density  of  8  x  10”  2  gm/cm2. 
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Fig.  4  (Continued)  —  Density  contours  as  in  Fig.  2  but  for  an  initially  uniform  4  x  10" 2  gm/cm 
-pinch  simulation  in  which  radiation  loss  terms  are  included.  6  nsec  and  peak  density  of 
.25  x  10"  2  gm/cm3.  (b)  14  nsec  and  peak  density  of  7.4  x  10"  2  gm/cm  .  (c)^16  liSec  l*8* 
density  of  1.7  x  10_1  gm/cm3.  (d)  18  nsec  and  peak  density  of  8  x  10"  2  gm/cm  . 


X  fO 


Fig.  4  (Continued)  —  Density  contours  as  in  Fig.  2  but  for  an  initially  uniform  4  x  1CT  2  gm/cm3 
-pinch  simulation  in  which  radiation  loss  terms  are  included.  8  nsec  and  peak  density  of 
.25  x  10"  2  gm/cm3.  (b)  14  nsec  and  peak  density  of  7.4  x  10“ 2  gm/cm3.  (c)  16  nsec  and  peak 
density  of  1.7  x  10“ 1  gm/cm3.  (d)  18  nsec  and  peak  density  of  8  x  10" 2  gm/cm3. 
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Fig.  5  (Continued)  —  Logarithmic  radiation  power  (ergs/sec-cm5)  contours  for  the  times  shown  in 
Fig.  4.  Each  plot  ranges  over  8  contours  with  3  contours  per  decade.  The  highest  possible  value  in 
(a)  is  4.6  x  1021  ergs /sec-  cm3  and  4.6  x  1022  ergs/sec- cms  in  (b),  (c),  and  (d). 


Fig.  5  (Continued)  —  Logarithmic  radiation  power  (erga/feec-cnr)  contours  for  the  times  shown  in 
Fig.  4.  Each  plot  ranges  over  8  contours  with  3  contours  per  decade.  The  highest  possible  value  in 
(a)  is  4.6  x  1021  ergs/sec—  cm5  and  4.6  x  1022  ergs/sec- cm3  in  (b),  (c),  and  (d). 
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Fig.  6  —  Total  line  radiation  as  a  function  of  initial  density. 
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